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Abstract 



o 
o 

£N| ■ The Volume Operator plays a crucial role in the definition of the quantum dynamics of Loop Quantum 

Gravity (LQG). Efficient calculations for dynamical problems of LQG can therefore be performed only if one 
has sufficient control over the volume spectrum. While closed formulas for the matrix elements are currently 
available in the literature, these are complicated polynomials in 6j symbols which in turn are given in terms of 
Racah's formula which is too complicated in order to perform even numerical calculations for the semiclassically 
important regime of large spins. Hence, so far not even numerically the spectrum could be accessed. 

In this article we demonstrate that by means of the Elliot - Biedenharn identity one can get rid of all the 6j 
symbols for any valence of the gauge invariant vertex, thus immensely reducing the computational effort. We 
use the resulting compact formula to study numerically the spectrum of the gauge invariant 4 - vertex. 

The techniques derived in this paper could be of use also for the analysis of spin - spin interaction Hamil- 
( ^^' tonians of many - particle problems in atomic and nuclear physics. 
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1 Introduction 



The volume operator ^ |5j plays a pivotal role in the definition of the quantum dynamics |1J |S| of Loop 
Quantum Gravity (LQG) [3]. Since the success of LQG depends on whether the quantum dynamics repro- 
duces classical General Relativity (GR) coupled to quantum matter in the semiclassical regime it is of outmost 
importance to know as much as possible about the spectrum of the volume operator. 

The volume operator has been studied to some extent in the literature |?1 181 IH1 HO] and it is well known that 
its spectrum is entirely discrete. However, so far only a closed formula for its matrix elements has been found. 
Unfortunately, not only is the formula for the matrix elements a complicated polynomial in 6j symbols involving 
extended sums over intertwiners, in addition the 6j symbols themselves are not easy to compute. Namely, the 
only known closed expression for the 6j symbols is Racah's famous formula which in turn involves fractions of 
factorials of large numbers and sums whose range depends in a complicated way on the entries of the 6j symbol. 
Accordingly, even powerful computer programs such as Mathematica or Maple run very fast out of memory even 
for moderate values of the spin labels on the edges adjacent to vertex in question. For instance, the current 
authors were not able to go beyond j = 3 when numerically computing the eigenvalues for a gauge invariant, 
four valent vertex, just using the matrix element formulas available in the literature. Thus, in order to make 
progress, analytical work is mandatory. 

In this paper, which is based on the diploma thesis we simplify the matrix element formula as given 
in |10) tremendously: Using an identity due to Elliot and Biedenharn we are able to get rid of all the sums 
over intertwiners and all the 6j symbols in the final formula, no matter how large the valence of the vertex is. 
The closed expression we obtain is a harmless polynomial of simple roots of fractional expressions in the spins 
and intertwiners, without factorials, that label the spin network functions in question. We reproduce the closed 
expression for the gauge invariant four - vertex which has been discovered first by de Pietri 

This formula should be of interest for a wide range of applications. First of all, it opens access to the nu- 
merical analysis of dynamical questions in canonical LQG. In particular, there is now work in progress aiming 
at extending the spectacular results of |12) from the cosmological minisuperspace truncation to the full theory. 
Possible first applications are alluded to in the conclusion section. Next, the techniques presented here could be 
of use for numerical investigations of convergence issues of spin foam models, see e.g. and references therein. 
Furthermore, our methods reveal that time has come to put LQG calculations on a supercomputer. Finally, it 
is conceivable that our formalism is of some use in the physics of many particle spin - spin interactions as e.g. 
in atomic or nuclear physics. 

The present paper is organized as follows: 

In section two we review the definition of the volume operator as derived in |2] and the closed expression 
for its matrix elements established in |10|. Knowledge of LQG £P is not at all necessary for the purpose of this 
paper which can be read also as a paper on the spectral analysis of a specific interaction Hamiltonian for a large 
spin system. 

The main result of this paper is contained in the third section where we derive the simplification of the 
matrix elements. At the danger of boring the reader we display all the intermediate steps. We do this because 
we feel that without these steps the proof, which in part is a complicated book keeping problem, cannot be be 
understood. The compact final formula is 147H . 

In section four we use our formalism in order to study the gauge invariant four vertex. The simplification 
of the matrix element formula now enables us to diagonalize the volume operator in a couple of hours for 
spin occupations of up to a 2j max « 10 2 . More efficient programming and compiler - based programming 
languages such as Lisp should be able to go significantly higher. Among the "spectroscopy experiments" we 
performed are the investigation of the computational effort, the possible existence of a volume gap (smallest 
non - zero eigenvalue), the spectral density distribution and the relative number of degenerate (zero volume) 
configurations. Among the surprises we find numerical evidence for a universal density distribution in terms of 
properly rescaled quantities valid at large spin. Next, there is numerical evidence for the existence of a volume 
gap at least for the four - valent vertex. Finally, it seems that the eigenvalues form distinguishable series just 
like for the hydrogen atom which provides a numerical criterion for the question which part of the spectrum 
remains unaffected when removing the finite size "cut - off" jmax ■ 

In section five we summarize our results and in the appendices we provide combinatorical and analytical 
background information which make the paper hopefully self - contained. 
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2 Revision of Known Results 



This section summarizes the definition of the volume operator of LQG and reviews the matrix element formula 
proved in |1U|. Readers not familiar with LQG can view the volume operator as a specific spin - spin interaction 
Hamiltonian for a many particle system. After some introductory remarks for the benefit of the reader with an 
LQG background we will switch to a corresponding angular momentum description immediately which makes 
knowlegde of LQG unnecessary for the purposes of this paper. 



In LQG typical states are cylindrical functions / 7 which are labelled by graphs 7. The graph itself can be 
thought of a collection £(7) of its oriented edges e which intersect in their endpoints which we call the vertices 
of 7. The set of vertices will be denoted by V(7). The cylindrical functions / 7 depend on SU(2) matrices h e 
which have the physical interpretation of holonomies of an SU(2)— connection along the edges e € E(-j). 

In 012] the operator describing the volume of a spatial region R, namely the Volume operator V(R)~f acting 
on the cylindrical functions over a graph 7 was derived as: 



(1) 



where 



V v . 



\ 



3! • 8 



2J e(e I ,e J ,e K ) qijK 



qijK = eijkXjX 3 jXfc 



(2) 



(3) 



(4) 



The sum has to be taken over all vertices v G V(j) of the graph 7 and at each vertex v over all possible triples 
(ei, e.j, etc) of edges of the graph 7 adjacent to v. Here e(ei, ej, &k) is the sign of the cross product of the three 
tangent vectors of the edges (ej, ej, e^) at the vertex v and we have assumed without loss of generality that all 
edges are outgoing from v. 

The X} are the right invariant vectorfields on SU(2) acting on the holonomy entries of the cylindrical 
functions. They satisfy the commutation relation [_Xj,-Xj] = —2 Su e lih Xj. The self-adjoint right invariant 
vector fields Yj := §Xj fulfilling [Yj, Yj] = i 6u e ljk Y r k are equivalent to the action of angular momentum 
operators J}. It is this algebraic property which we use in order to derive the spectral properties of the volume 
operator: It turns out that the Hilbert space of LQG reduces on cylindrical functions over a graph 7 to that of 
an abstract spin system familiar from the theory of angular momentum in quantum mechanics. There are as 
many degrees of freedom as there are edges in 7 and furthermore we can diagonalize all the V v ,-y simultaneously 
as they are obviously mutually commuting. Hence, in what follows familiarity with LQG is not at all necessary, 
abstractly we are just dealing with an interaction Hamiltonian in a many particle spin system. 

We can therefore replace: 



qijK = 



jt jj jk 



(5) 



Using furthermore the antisymmetry of tijk and the fact that [J}, Jj] = whenever I 7^ J we can restrict the 
summation in @ to I < J < K if we simultaneously write a factor 3! in front of the sum. The result is: 



KJ<K 



(6) 
(7) 



Now the following identity holds: 

Uik JIJ'jJk = \ [{Ju)\{Jjk) 2 ] (8) 

where Ju — Ji + Jj. This relation can be derived by writing down every commutator as |"( Jjj) 2 , ( Jjk ) 2 j = 
S |^( J) + J}) 2 , (Jj + J^) 2 j , using the identity |^a, fecj = ^a, fej c + fe|^a, cj for the commutator, using the angular 
momentum commutation relations 11041 and the fact that [jj, Jjj = whenever I ^ J. 
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We may summarize: 



K,7= £ e ( e i, e J> e i<) qijK\ (9) 

y K.KK 

(10) 

where q IJK ■= [(Ju) 2 , (Jjk) 2 ^ and Z = J. 

Unless announced differently we will study the operator 

QIJK := [(Jij) 2 ,(Jjk) 2 ] (11) 

in the following. 

2.1 Matrix Elements in Terms of 3rv/-Symbols 

Now we can apply the recoupling theory of n angular momenta to represent qijK in a recoupling scheme basis 
using the definitions 1A.1L 1A.2L 1A.3L given in the appendix. 

We will do this with respect to the standard basis l|A.20 . where we can now easily restrict our calculations 
to gauge invariant spin network states, by demanding the total angular momentum j and the total magnetic 
quantum number M to vanish, that means we will take into account only recoupling schemes, coupling the 
outgoing spins at the vertex v to resulting angular momentum 0: 
In terms of the recoupling schemes these states are given by: 

\g(IJ)jj = 0M = 0>:=\g(IJ)> (12) 

where we introduced an abbreviation, since the quantum numbers j j = M = are the same for every gauge 
invariant spin network state with respect to one vertex v. 

We will now represent qua '■= ^{Ju) 2 , (Jjk) 2 ^ in the standard recoupling scheme basis of definition l|A.2fl 

where \a >:= |a(12) >, \a! >:= \a'(12) >. 

The point is that by construction a recoupling scheme basis \g(IJ) > diagonalizes the operator 
(G 2 ) 2 = {Juf = {Ji + Jjf that is 

(G 2 ) 2 |g(/J) >= g 2 (U) ( 52 (/J) + l) \g(U) > (13) 

Furthermore every recoupling scheme \g(IJ) > can be expanded in terms of the standard basis via its expansion 
coefficients, the 3nj-symbols given by definition I A . 31 in the appendix. So it is possible to express 

< a\quK \a > = < a\ [(Ju) 2 , (Jja-) 2 ] \a > 

= < aKJufiJjK) 2 ]^' >-< 3\{JjK?{Juf\a > 

= J2 MU)(g2(IJ) + 1)[< a\g(IJ) >< gu\jj K \a > - < a\j'j K \gu >< g{IJ)\a >] 
S(i J) 

E g2(U)(g2(IJ) + l)g2(JK)(g 2 (JK) + l)<g(IJ)\g"xg(JK)\g">x 

g(IJ),g(JK),g" (12) 

x[< g{IJ)\a>< g{JK)\a > - <a(JK)\g>< g(IJ)\a >] 



= E 



s"(12) 



£ ,g 2 (/J)( 5 2(/J) + l) < g{U)\g" >< g(U) 



a > x 



x £ g 2 (JK)(g 2 (JK) + 1) < g(JK)\g" >< g(JK)\a' > 

3(JK) 

-[*<=> a'} (14) 

which is again an antisymmetric matrix possessing purely imaginary eigenvalues (we could alternatively consider 
the hermitian version by multiplying all matrix elements by the imaginary unit i). Here we have inserted suitable 
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recoupling schemes \g{IJ) >, \g(JK) > diagonalizing (Ju) 2 and (Jjk) 2 and their expansion in terms of the 
standard basis |<?(12) > by using the completeness of the recoupling schemes \g (IJ) > for arbitrary I 7^ J 1 : 



So we have as a first step expressed the matrix elements of qijK in terms of 3nj-symbols. 

2.2 Closed Expression for the 3nj-Symbols 

The 3nj-symbols occurring in 11 H can be expressed in terms of the individual recouplings implicit in their 
definition. 

2.2.1 Preparations 

In |10| the two following lemmas are derived: 
Lemma 2.1 Contraction on Identical Coupling Order 

< g(U)\g > 

=< 92(31, Jj), 93(92, ji), ..,gi+i(gi,ji-i),gi+2(gi+i,ji+i), ■ ■,gj(gj~i,jj-i)\ 

\92(jl,j2), g'3(92, 33), gi+i(9l , 3l-i), 9l+2(9l+l, Jl+i), ■■■■,9j(9j~iJj) > Ggj.g'J-Sgn-i.g''^ 

Lemma 2.2 Interchange of Coupling Order 

< 92(31,32), ■■,9'k(9'k-i,3k),9'k+i(9'k,3k+i),9'k+2(9'k+i,3k+2)\ 
\g2(ji,32),-;9K(gK-i,3K),gK+i (gx , 3K+2) ,gK+2 (gx+ 1 , ji<+ 1 ) > 

= < 9K+l(9K , 3K+l), g'K+2(g'K+l, 3K+2)\9K+l(gK, 3K+2), 9K+2(9K+1, 3K+l) > ^g' 2 g 2 ^g' 3 g3—^g' K gK ' 3g' K+2 g K +2- 

2.2.2 Closed Expresion for the 3nj-Symbols 

Now we can reduce out the 3nj-symbol. In what follows we will not explicitly write down the ^-expressions 
occurring by using lemma I?TT1 and lemma [5^] but keep them in mind. 

Collecting all the terms mentioned in |10| one obtains the following equation for the 3nj-symbols: 



1 The summation has to be extended over all possible intermediate recoupling steps g2, ■ ■ ■ ,g n -i that is \j r —jq\ < 9k(jq,jr) < jq 
allowed by theorem IA.1I given in the appendix. 




(15) 



< g{U)\g(l2) > 



^2 < 92(ji, jj), 53(32, ji)\h2(ji,ji),g 3 (h2,jj) > 
< 9z(h2, 3j), 94(93,j 2 )\h 3 (h 2 ,j 2 ), g4(h 3 , jj) > 



^2 < gj-i(hi-2,jj),gi(gi-i,ji-2)\hi-i(hi-2,ji 3 ),gi(hj-i,jj) > 

h I-l( h I-2,3l-2) 

X < g I (h I - 1 ,jj),g I+1 (g I ,j I - 1 )\g' I (h I - 1 ,j I - 1 ),g I+1 (g' I ,jj) > 

x < gi+i(g'i,jj),gi+2(gi+i,ji+i)\g'i + i(g'i,ji+i),gi+2(gi + i,jj) > 

x < gi+2(g'i+i,jj),gi+3(gi+2,ji+2)\g'i + 2(gi+i,ji+2),gi+s(gi+2,jj) > 

x < gj-i(gj-2,jj),gj(gj-i,3J-i)\g'j-i(g'j-2,jj-i),gj(g'j~i,jj) > 

x < h 2 (ji,ji),h 3 (h 2 ,j 2 )\g 2 (ji,j 2 ),h 3 (g 2 ,ji) > 

x < hz(g' 2 ,3i),hA(hz,33)\g'z(g' 2 ,3z),hi(g' z ,3i) > 

X < h 4 (g 3 ,j I ),h 5 (h4,j4)\g'4(g'3,34),h5(g'4,jl) > 

x < hi-i(g' I _2,ji),g'i(hi-i,ji-i)\g' I _ 1 (g' I _ 2 ,ji-i),g' I (g' I _ 1 ,ji) > 



(16) 



3 Simplification of the Matrix Elements 

3.1 3rr/-Symbols Expressed in Terms of 6j-Symbols 

It is now obvious that we can express l|16jl via the 6j-symbols defined as in l|119|l : 

< 312(31, 32), 3(312, 33)\h3(j2, 33), j(ji, J23) >= 

= [(2j 12 + l)(2j 23 + 1)] 1 (_i^+i a +i a +j I 3i 32 312 1 

L J3 J J23 J 

Here we have used the fact that j\ + j 2 + j 3 + j is integer. Now the definition of the 6j-symbol in terms 
of Clebsh-Gordon-coefficients (CGC) come into play. Because of the properties of the CGC we can change the 
order of coupling in every recoupling scheme in l|16|l taking care of the minus signs we create: 

< 3 12 O'l ,32), j (jl2, 33 ) b'23 (32 ,33), j (jl , J23 ) > = 
= < j 12 (j 2 , Ji), j(j 12 , j 3 )\j 23 (j 2 , 33), j(3l, 323) > 

= ( _ 1) i 12 - J - 1 -i 2( _ 1) i 23 -i 2 -ia <jl2 (j 2 J 1 )J{j 12 ,j 3 )\j 23 (J 3 ,j 2 ),j(j 1 ,j 23 ) > 

= ( _ 1) «2-«-,- 2 (_ 1) J 2 3- 3 -2- J -3(_ 1 y-J I2 - J -3 <j 12 (J 2 J l ),j(j 3 ,j l2 )\j 23 (j 3 ,j 2 ),j{j 1 ,j 23 )> 

= ( _ 1)J l 2 -,l-, 2( _ ir 3-, 2 -« ( _ 1) ,-, 12 -,3 ( _ 1) ,-,l~ J2 3 < j 12 (j 2 , ^0,^03, il2)|j23(j3, j 2 ), J(i23, jl) > 

In this way we are able to change the coupling-order in 1161 to get the order required for a translation into the 
6j'-symbols. With these preparations we are now able to express l|16> in terms of 6j-symbols: 
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< g(I.J)\g' (12) > = 

J2 (_l)-JJ-JJ+92(_l)h2+J.7-S3(_l)W+JJ+Jl+93 ^(232 + l)(2/l 2 + 1) <{ V ^ ^ 



J, 2 



X 

h-3 



J2 (_l)-Jj-''2+93(_l)''3+JJ-94(_l)JJ + '>2+J2+S4 v /(2 S 3 + l)(2/l 3 + 1) | V ^ ^ 



X 



^ (_l)-JJ-h/-2+9/-l(_l)h/-l+JJ-9J (_l)JJ+hj-2+J/-2+9f % /(2ff/_l + l)(2/lj_l + 1) | 

(_l)-ij-hi-i+9j(_l)jj+9i-9i+i(_l)ij+fcj-i+ji-i+9i+i ^(237 + l)(2pj +1) I 
(_l)-jj-9j+ff/+i (_i)JJ+9j+i-9i+2(_i)ij+Bi+i/+i+9/+2 s J{2gi +1 + l)(2g' I+1 + 1) j 
(_l)-JJ-9i + i+9J+2(_ 1 )Jj+9i + 2-9j+3(_ 1 yj+9r + i+Jr+2+9j+ 3 y(2 g/+2 + l)(2g^ +2 + 1) I 



jj Si-i 



91+1 Sj 
JJ 9/ 91+1 

ji+i 9i+2 g' I+1 



03 Sj+1 9/+2 
ji+2 3/+3 9/ +2 



3J 9j- 2 9J-1 
9J Sj-i 



X (_l)-W-9j-2+9.7-l(_l)J.7+9j-l-9.7(_l)J.7+9j-2+JJ-l+9./ v /(2g J _ 1 + 1)(2 9 ' J _ 1 + 1) j ^ J j 

X (_l)i/+92-h3(_l)Jj+ii+J3+h3 s /j2M + l)(2g' 2 + 1) | ^ ^ 

X (_l)-9i-ij+fc3(_l)«+9 3 -'^(-l)3J+93+i 8 +h4 A /(2/ l3 + 1)(2^ + 1) j ^/ ^ ^ 3 

X (_l)-93-J/+h4(_l)JJ+9i-h5(_l)J/+93+J4 + h5 v /(2/ l4 + l)( 2 g£ + 1) j ^ ^ ^4 

x [- 1 )-3i-2-Ji+hi-i ( _ 1 yi+g' I - 1 -3i { _ 1 yi+a' I - 2 +ji-i+s' I ^ 2hj _ 1 + i)( 2ff ^_ 1 + 1) | ^ 9 ^2 ^-1 

(18) 

This is the complete expression of Q16|l with all the exponents written in detail which are caused by the reordering 
of the coupling-schemes while bringing them into a form suiteable for 111911 . We want to emphasize that we 
have the freedom to invert the signs in each of the exponents of l|18fl when convenient for our calculation. 

3.2 The Matrix Elements in Terms of 6j-Symbols 

Taking a closer look at 1141 , a basic structure contained in the matrix elements of the volume operator appears: 
E 92(U) ( 92{U) + 1 ) < g(IJ)\g" (12) >< g(IJ)\S(12) > (19) 

Using GBJ) we now express the 3nj-symbols occurring in Ijl9^ via 6j-symbols. For < g(IJ)\g" (12) > we use 
h,2 . . .hi-i as intermediate summation variables and for its (—1) -exponents the sign convention we chose in 
11811 . For < g(IJ)\S(12) > we use . . . fcj-i as intermediate summation variables and for its (—1) -exponents 
the negative of every exponent in l|18fl . since every exponent is an integer number. Writing down carefully all 
these expressions most of the exponents can be cancelled. 
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The result of this is (using the abbreviation A(x,y) = y (2x + l)(2y + 1) ): 
E 92( /J ) ( 92(U) + 1 ) < g(IJ)\g"(12) X g(IJ)\d{\2) > = 



9(1 J) 



E 92 {I J) (p 2 (JJ) + l 



9(U) 



X 



X 

h 3 



X 



E (-i)*ag»,*o{£ J £} E (-i)-M( ff2>fe) { j £ g 



( _ l) 9'UiA( gi+1 ,g''){ jj 9 " (-l)-^ +1 A( Sl+lj a /+1 )( ai 91+1 

v I \m+ua I+ u | gi+2 g ii +i J v ; w+i' \ 9i - +2 a /+ i 



(-l)fcs+tfA(hs,i#){2 g g j (-l)-*s-«»A(fc3,a9) { 

( -l)^'^ 4lff «){jgJ|} 



(-l)-*4-<*4yl(fc 4)04 ) 



jl 02 &3 

J3 ^4 «3 

JI a.3 ki 

J4 ^5 «4 



x (-i)V-i+rf_ lA(h v )( J'/ S "-2^ X 1 (-l)- k '-i-MA(k Wl , w )|/' a/_ a fcj_! 



(20) 

Up to now we have only made a translation between different notations. The reason for doing such an amount 
of writing will become clear soon: Using identities between the 6j-symbols it is possible to derive a much shorter 
closed expression for the matrix elements by evaluating step by step all the summations in 120H . 

3.3 A Useful Identity 

Before we can start the evaluation we want to derive an identity, which will be essential. We want to evaluate 
the following sum: 

fWn) - Vf2„, + 1) + 1) { * f f } { * f f } (21) 

^ I 33 H 323 ) [33 3i 323 J 
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There exist numerous closed expressions for 6j-symbols, whose entries have special relations. They are much 
more manageable than the general expression given in the appendix I12UH (Racah formula). One of them reads 
as HU, p.130: 

f a b C 1 = 1)a+i; + c +i 2[6(6 + l)+c( C +l)-a(q+l)] (22) 
l lc6 J [26(26+ l)(26 + 2)2c(2c + l)(2c + 2)]3 

Using the shorthand X(b, c) = 26(26 + 1)(26 + 2)2c(2c + l)(2c + 2) and the fact that a + b + c is integer, we can 
rewrite the equation to obtain 



a(a + 1) = (-!)<•+"+* | ° ^ } X(6, c)* + [(6 + 1) + c(c + 1)] 
Putting a = J23 and inserting 1231 into 12 1 1 one finds for F(jii,j'i2) for any 6, c : 



(23) 



323 

6(&+ i) + c ( c+ i) f * j2 il2 1 f * J2 & 1 

(2J12 + 1) ^ V I J3 34 J23 J Us 34 J23 J 



where we have used the orthogonality relation 11241 for the 6j-symbols. Let us take a closer look at the three 
6j-symbols of I on the right hand side of 1)24(1 . We now apply some permutations to the rows and columns of 
the first 6j-symbol within / which leave this 6j-symbol invariant, see (|122( . I)123( , After that we have for /: 

i = E(- 1 ) J23 (2j» + 1) { 6 6 1 W h 32 J - 2 } { 31 J - 3 ' 12 

4^ I C C J 23 J [ 33 34 323 J [ 33 J4 323 

J23 



ji jl 1 1 f ji ia jl2 j j jl j2 jl2 
34 34 323 J 1 J3 jl j23 J 1 33 jl j23 



jl2 


jl 




{ js 


jl2 


312\ 


jl 


3 12 






34 


j* J 



Now we fix 6 = ji , c = ji and can evaluate: 

1 = E(- 1 )' 23 (2j23 + l)' 
J23 

_ /_^\-0l2+Jl+Jl+J3+J2 

Here we have used the Elliot-Biedenharn-identity 1126H . Inserting this back into 1241 yields: 

F(ji2,j' 12 ) = I(-l) J1+ ^(-l)-« 12+jl+ ^^ 

, Jl(jl+l)+j4(j4 + l) ) - 

(2j 12 + 1) J12 ^ 

Using that, by definition of the 6j-symbols l|119[l . (J\ 2 + ji + j 2 ) and (js + 1 + j 4 + j' 12 ) arc integer numbers, we 
can invert their common sign in the exponent of the (—1) . After summing up all the terms in the exponents 
and performing some permutations on the arguments of the 6j'-symbols according to (11228 . 1123H we obtain the 
final result for F(ji 2 , ji 2 ): 



f jl2 


jl 




u 


312 


fl2\ 


I jl 


jl2 






34 


j* J 



F{ji2,j'i 2 ) := J2^ 23 + X ) ■ ?23 ^ 23 + X ) 

323 

1 / J \l'i4-i'o-l-1': 





32 


jl2 1 




72 


j; 2 i 


I 73 


34 


j23 J 




74 


j23 J 





jl 


*} 




.74 


jl2 1 




ji 2 






.712 


j* J 



Jl(jl + l)+j4(j4 + l) 
(2jl2 + 1) 



+ wi ,:rr x; w u (25) 



with X(ji, j 4 ) = 2ji(2ji + l)(2ji + 2)2j 4 (2j 4 + l)(2j 4 + 2). 
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Remark By the integer / positivity requirements of the factorials occurring in the definition of the 6j-symbols 
of see (Hill we can read off restrictions for ji2,j'i2, namely the selection rules 

f Jl2 ~~ 1 

j'12 = < 312 

{ Jl2 + 1 

3.4 Precalculation 

After these preparations we can now go into 'medias res': We will carry out all the summations in (121)1 . Step by 
step we will only write down (with reordered prefactors) the terms containing the actual summation variable, 
suppressing all the other terms and sums in I20H . 

We start with the summation over g(U) (using again the shorthand A(x, y) — \J (2x + l)(2y + f ) and the 
fact, that A(g2, /12) ■ A(g2, k 2 ) = A(h2, fe) ■ (2(/2 + 1)- Additionally we will frequently use the integer conditions 
E3): 

• Summation over g 2 : 

^.^(-D^E^+DPto+DiS* £}{£££} = 

92 



I122J123I 





31 


h 2 1 




h 


k 2 1 


I 93 


3J 


92 J 


I 93 


3J 


92 J 



EH 



A(h2,k 2 )(-l) h2 - 



-(-V)3l + 31 +h 2 +g 3 +jj+l+h2 X(ij a \\ [ 31 31 h 2 \ f 93 jj h 2 \ 
2 ( 1 *VI,3J) \ 1 k2 ji j \ 1 k 2 3 J J 



+ Sh * 

A(h2,k2)^(~ir jl - J1+Sj+1 X(j I ,jj)i l j l £ J 



( _ 1)9 , J 93 jj M + [ + 1} + !) ] 5ft2fc2 



M 2 



• Summation over 33: 
A(h 3 ,k s )(-l) h 3- k iJ2 



M 2 (-l) 5 



93 jj ^2 
1 k 2 jj 



(293 + 1) 



3 J h 2 93 1 f jj *=2 93 
J2 94/13/1 J2 94 ^3 



{ ^E(-i)«(^ + d { ? £ * } { j g } { % } 



™ a ^3,*3) ( -i)—3 { m 2 e(-d-( 2 93 +d{ £ £ i } { s 2 £ } { £ s 5 } 

:- 5 h2'=2E( 2fc 3 + l)(293 + l) 



93 
N 



(2fe 3 + l) 



(JJ 


94 


fa 1 


f JJ 


94 


5}} 


1 J2 


fa 


93 J 


1 J2 


fa 





Ufa. !•:;)(- I.)" ; 1 1 <| M 2 (-l)- 

iV 



l 12b l 124 l 4fh, ^V-1^3- fc 3 J Mo(--11-( h 3+i.7+i2+JJ+94 + l + fc3 + /l2 + fc2) / ft 3 JJ 94~l f 1 fa fc 3 

JJ fc 3 1 J 1 32 k 2 h 2 



'(2k 3 + l) 5h2h25h3k3 



II: JIIJJII j~1 



M2A {h3 , k3 )(-ir^ +k ^{ T t ll h k l}(-D 9 *{ 9 ?i J 3 $j}+*w*s*, 



M 3 

Note: that is in principle the same term with the same index-order as we got from the summation over 

92- 
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Generally we have for the summation over gi for 4 < i < I 



Ah . . . { 9 l l J _^ ^r 1 } + N8 h2k2 . . . 

jj hi-i 9i \ [ jj ki-i gi 



x (2 ff < + 1) 



ji-l 3*+l hi J 1 g j+ i fcj 



A(h j ,fc < )(-l)' l *- fc ' x 



M 2 ...M l _ 1 ^(-l)^(2 9l + l)| * £ ^j" 1 

+7V5 h2fc2 ...5 hi „ lfci _ 1 ^(2 9l + l)] JJ £ 

^7" ^ ji— 1 y-i+i n z J ^ Ji — l in+i K t 

A(h t ,ki)(-l) h >- k > x 

A/ 2 ... Mi_! J2 ( 2 ft + 1 ) 



Jj fti-i 9i \ f Jj fc i-i Si 
ji-1 Si+1 ft i J I 9i+l fc i 



Jj hj-i Si \ / JJ fci-i Si 



nzana 



j'j jj 1 1 j jj ^ Si+i \ r jj s,+i fcj 

ki-i Si J 1 Ji-i Si hi-i J 1 Ji-i fci-i Si 



(2fcj + 1) 



I Ji— l 



jj ^ Si+i i r jj ffi+i 

Si hi-j J 1 jr'i_i fei_l 



umnm 



integer integer 

x I Mo M- 1 f— 11 — J+Ji-i+ij + ffi+i + l + ki + hi-i + fcj_i) f hi jj Si+i \ ( 1 fei 
[ v ; \jj ki 1 J \ ji_i 

+N8 h2k2 . . . 8 hi _ lki l 8 h . ki 

M 2 ...M i _ 1 A(/ liJ fc i) (-i)- i *- i+fc *- i+& - i+i { : 'V ^ ^(-^{t £ £} 

. ' 

+NS h2k2 . . . &h i _ 1 k i _ 1 &h i k i 

• For the summation over gi+i 

we have the same terms as above with a slight difference in one index (underlined): 



Ili.-'IIJ-'llJl 



m 2 . . . Af/(-i)«+» { 9 y x £; J; } + • ■ • ^ 



jj fei Sj+i 



x (2si +1 + 1) ( h hl f +1 , , ■ 

\ Jl+l S/+2 "7+1 J 1 J/+1 Si+2 fcj+1 



M 2 . . . M 



9i+2 jj h I+ i 
ki+i jj 



Ml+l 

+N8 h2k2 . . . 8h I k I $h I+1 k I+1 
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• Summation over gi for I + 2<i<J — 1: 



A(/ ll ,fc l )(-l)' 1 *- fe -£ 



Af 2 ...Af/M/ + i...iQi_i(-l) 9 ' 
x (2 ffj + 1) 



l fci-i jj 



+ N5 h2k2 . . . S h ._ lk ._ 1 



ji g i+ i hi J 1 ji g i+l ki 



M 2 ... MjM I+1 . . . M,_i A(hi, ki){-l)-^+ h 



1 k, kiU f ( 1)S! + 



i / 9i+i jj hi 
1 h jj 



Here we keep the following notation in mind: 



hi 


= Ji 


fci 


= 31 


J <i < N : gt — g'( — en 




h 2 


= h 2 (ji,ji) 


k> 








h 3 


= h 3 (h 2 ,j 2 ) 


k 3 


= ka(k 2 ,j 2 ) 






hi-i 


= /l/-l(/l/-2,j/-2) 


ki-i 


= fc/_l(fc/_ 2 , j/-2) 






hi 


= g" 


ki 


= a/ 






hi+i 


a 

= 9i+i 


ki+i 


= 






hj-i 


= g'J-i 


kj-i 


= aj-i 
















(26) 



We have now carried out completely the summation ove r g(7 J) in (EDI- After this we write down the re- 
maining terms of this summation and the terms of 12U1 which have not taken part in the summation yet. 

So we end up with (do not get confused about the notation A(x, y) = U (2x + l)(2y + 1) whereas A, Ai are 
abbreviations form certain terms): 

52 92(32 + 1) < g(U)\g" >< g(U)\S > = 

9(U) 



part I 



E 



h 2--- h /-l 
^2-^1 — 1 



n=2 m=I+l 



n M n rr M m x (-1)^ \ 7 33 



part II 
J-l 

n 



X remaining terms 



E 



J - 1 . v-'^-l 
3 J— 31 — Z J „=i Jn-L m=J+ i 3n 



Ai<—>n=3,A2<—>-n='i,A3*—fn=5,. . . , Aj _2 *—n%=I 



x (-l) +1 A(/ l2 ,fc 2 ) 



ii Ji ^2 
1 jj 



n A(/ ln ,fc„)(-i) h "-i+ fe "-i+ i ^ n 1 

n— 3 
J-l 

]^[ A(^,a„)(-l) 

x (-!)-■' { y a ^ g " j .- x } + [ji(ji + 1) + jjuj + 1)] n 



jn—1 hn—i h n 



9^'_l+a„-l+l J J™ Sn-1 9n 

1 a n a n _i 
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B 3 



(-l)- a ^(fc 2 ,a 2 ) 



J2 Jl &2 
J2 &3 «2 



C 2 

(-1)-M(fc4,«0(^ I 3 M 
v ' v ' |_ 34 k 5 C14 J 



C 3 



( _ 1) hi- 1 -fc/-i+ S i'_ lj4( ^_ li 2^_ 1 ) 



ff"_2 hi- 
ll -l Sj ffj_ 



(_1)-»i-iA(Ai_i,o/_i) 



(27) 

After noting this intermediate result we finally have to execute the remaining summations of Ij27ft , namely the 
summations over h 2 , . . . , hi-i and fc 2 , ■ ■ ■ , fcj-i (leaving out the signs (— l) 9i _a ' , since they will be cancelled, 
as we will see, due to the occurrence of S g ii a . -terms in the following calculations). 



First we do this summation for the 

• first step 

summation over /i 2 : 

A 



part I of 



h 2 



31 J2 9 2 



f fc 3 


kz 


32 \ 


\3! 


k 3 


f} 


I 31 


92 


31 J 


\ h 3 


3i 





1 J 1 ji 2? 2 ' JJ J I As 3i 1 J 



(_l)tf+i+*3-*-»-*s A (A3,*a) { j J f } A(fc 2>ff2 ') { * g * } 



El 



summation over fc 2 : 



Oi 



^A (fc2 , 9 »){^ k ;t, £ \ 



k 2 



31 92 3/ 



31 3i k 2 

32 k 3 a 2 



E(2 fc2+ l)( 2a2+ l) 
(2a 2 + 1) t 



fir 


32 




fir 


32 


a2 1 


I fe 3 




fc 2 J 


1 fc 3 


.7/ 


fc 2 J 



EH 



9 2 "2 



(28) 



S 2 »2 



(29) 
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second step 

summation over h 3 : 



A-2 



(_l)92-Ji-i2+fc3 J 4( fc3i9 »)A(/ l4 ,fc 4 ) ^(-l) h 3(2A 3 + i) 

>>3 



1 n 


fe 3 


f} 


f i3 /13 


h 4 \ 


/ h 9$ 


h 3 1 


\ h 3 


J/ 




1 1 fc 4 


*; 3 J 


\ 33 h 4 


S3' J 



ll---1IJ.il 



E*. (-!)*»(«. + l){f J £}{ 



&3 &4 J3 1 f 3 2 ' S3 
?»4 ?»3 1 J I /l4 j/ ^3 



(_l)«s 



+ l + *4-il -32 -33 - 



EMt^-l^+j! + 1+93 +k3+9%+h4+u) J fc 4 *=3 is 1 f 3l k 4 g'i \ 

\ a'i g'{ 31 J I h 4 1 J 

DS1 ( — l) fc 4 +J7 +1 + 3;" - fc 3 -9 2 ' - ^4 ~33 ~3I 

^ 4 , fc4) {^J; 4 f }A (fc3 , 9 3'){ fc4 fc3i3 



&4 &3 j3 1 f 31 ki g 3 
92 S3 J* J I ^4 1 



9 2 ' S3' Jl 



D 2 



E-2 



(30) 



summation over £3: 

E-2 



Co 



E A ( fc 3-S3) 

fe 3 



i/ a-2 k 3 
33 ki a 3 



ki k 3 33 ] „/, \ 

(2a 3 + 1) 



9 3 03 



E(2fc 3 + l)(2a 3 + l) 



^3 



/ «2 


J3 


9'n 


/ 0.2 


33 


o 3 \ 


\ fc 4 


31 


k 3 J 


\ k 4 


31 


k 3 J 



S3 "3 



(31) 



• third step 

summation over /i 4 : 



D 2 



B 3 



j2 (-i) 9 3 

/14 



- h -A(h 4 , ki) { « J f } (-IJ^+Ufa, fc 5 ) { * 



J4 ^4 

fcs fc 4 



— ^4„4'){-g^ 



(-l)S3-«-J2-J 3 +fe4 A ( fc4i g») A ( ft5)fe5 ) (2^4 + 1) 

^4 





/i 4 




f j 4 h 4 h 5 \ 




S 3 ' 


h 4 1 


\ ^4 


J/ 


¥} 


1 1 k 5 k 4 J 


I J4 


/IS 


s 4 'i 



112211231 



{? 


fc 4 




f fc 4 


k 5 34 \ 


/s 3 ' 


,/4 


sn 






hi J 


\ hs 


hi 1 J 


\ ^5 


.// 


ft4 J 



IMl ( _ 1) -(fc 5 +i / + l+g 4 '+fc4 + 93' + h5+i4+i/) J fc 5 fc4 j 4 l f jj % S4M 

\ 9 3 9i 31 J I h 5 ^ ji 1 J 

I(_l)fc5+i/+l + 94 _ ' c 4-g3'-'l5-i4-jj . 



k$ ki ji \ J fc 5 g 4 ' 
S3 s" J I ^5 j/ 1 



(-1)94 



+ l + fe 5 -J! -J 2 - 33 -34. - 



h5 A(h 5 ,k 5 ) 



31 k 5 g'l 
h 5 31 1 



E3 



(32) 
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summation over ki\ 

E 3 C 3 



&4 



(2a 4 + l)£ 



/ «3 


34 




/ 0.3 


34 


a 4 "1 


1 fc 5 


ii 


fc 4 / 


\ k 5 


31 





In this way we successively carry out all the summations until the last step: 

• last step 

summation over hi-i: 

D1-3 

y (-i)s"-i+ 1 + fc x-i-E»=i^-''x-M(ft / _ 1 ,fe I _ 1 )/ h g'j' 2 1 x 



(33) 



(-l)«i'-»-Si=?i»+*r-i i i(fc I _ 1 ,fl?_ 1 )i4(fc / ,fc / ) x 

x E + { *, *-> *- } { v V £, } { if:, 1? } 

UM ( _ u -(ki +ji + l + gi_ 1 +k I _ 1 + g > I '_ 2 +h I +j I _ 1+jl ) ( ki if-iHj/ kj g'{_ x 

lMl(_i\ k i+ji + 1 +3i-i- k '--i--3'i'-2- h i-n-i-ji / fe J" 31-1 \ / W fc / g'l-x 



Dj-2 Bf-2 



(34) 



summation over 

Bl-2 Cj_ 2 



fe/_i 



jl 0-1-2 1 _ 



™a Afa - 1 ' 9 "- l) T (2k^ + i)(2 OI _ 1 + 1) ( a '- 2 ^r 1 f-i ) { a ^ 'V 1 

(2a 7 _i + l) ^ ; \ a/ J/ fej-i J \ ai 3i 

fej — 1 



CP*,, 



5 Sx'-i a i-i 



(35) 

Here we have used that <?"_ 2 = ai-2 resulting from the summation over ki-2 and additionally the fact 
that ai = ki from l|2fi[l . 



1G 



At the end of the part I-summation over h 2 • ■ • hi-i,k 2 • • ■ fcj-i we can now summarize the remaining terms 
of 1271 resulting in the term D1-2 from the summation over hi-i in 1341 1 

With <?"_i = aj-i from 13511 , hi — g" , ki — ai from Il2(it we obtain the drastically simplified formula: 



E 



n=2 m=I+l 



n M n n x (-1)^ <• 7 jj 



x remaining terms 



(-i)«'-i+Vir- ff 'Vir E ™ =i ^(^ai) 3* 



ji a/ ai-i 



M m x (-1)- ^ ^ 



m = J + l 



a J JJ 9j-i 
JJ 



(36) 



Now the summation for 



part II in 1271 is the last task in order to complete our calculation: Let us write down 



this expression (again suppressing the ( — l) 9i "'-signs) with the product 11^=2 ^h n k n evaluated (this cancels 
the summation over k 2 and all the exponents in (— l) hi ~ ki ): 



E 



&2 ■ ■ kj_ 1 



n n s ^ n 



X remaining terms ■ 



[Mh + 1) + 3 Ah + 1)] £ x M^aS) { % l l 3 J?, } 

h 2 ...h I _ 1 s ^ , 

y ' yi, \j3 h 4 g'l 



A{h 2 ,a 2 ) 



ji 31 h 2 
h kz a 2 



Oi 



A(h 3 ,a 3 ) 



3l «2 ^3 
33 ki a 3 



C 2 

A(h 4 ,a 4 )\ j ' I 3 h * 
1 34 k 5 0,4 



y 1 11 ' 31-1 9i 9i-x 



y4(ft/_i,a/_i) 



3 1 »/-2 ^1-1 
or a-l-l 



(37) 

Looking at 1)3 1\ we can see, that every summation gives rise to an orthogonality relation between the 6j- 
symbols as foll ows : 

Every sum in 137^ has the following form (again using the conventions in \Y2(il ): 



l 3n n n+1 g n 



> A(h n , On) 



31 ttn-l h„ 

jn hn-\-X 



A ^^i) y [2h 1)( 2 +1) \ji <_i h 



31 a„_i h n 

jn 



EH 



One has to start with the summation over h 2 . This gives 5 g q ao . 
Secondly the summation over fa 3 is carried out by using the result S g /_r a2 
way one can step by step sum over all hi up to hi-i. 
The final result for part II in l|27jl is: 



(38) 



from the summation before. In this 



E 



N\[5 h 



J —1 

remaining terms = + 1) + + 1)1 JJ 



(39) 
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Now we have solved the problem posed in 11201 and can write down the remarkably simplified expres- 
sion by using the results of (13611 . 1391 . Note, that we have inserted ^-expressions coming from lemma 



£ 92(H) ( g 2 (IJ) + l)< g(U)\g" (12) X g(IJ)\a(12) > = 



9(U) 



J-l 

n V(2<)(2a„ . + !)(-! 



J— 1 , (a a 

9„-l+<»«-l+l ] J™ 5n 



X 

n=I + l 

1-1 



1 a n a n _i 

J J n=2 n=J 

+ 1) + jj(jJ + !)] II ( 4 °) 



For configurations (/, J) one has to take all terms of I40H . which are in suitable limits, e.g. if J = I + 1 then 
j-i 

the product J] ... is not to be taken into account. Notice that for special configurations I < J certain terms 

n = J + l 

J-l 

drop out, e.g. if J = / + 1 then } J is not taken into account. 

n=I+l 
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Let us for clarity explicitly discuss four special cases of the edge-labelling I, J, namely |l = 1, J arbitrary 

I = 1, J = 2] , \l = 2, J arbitrary] and [ I = 2, J = 3] . 
We will display for every case the parts remaining from )t . 
Note, that again A(x,y) = y' (2a; + l)(2j/ + 1), and the conventions of i'2(>l are kept in mind : 



3.4.1 



I = 1, J arbitrary 



E ( fla(* -0 + 1 ) < 9(U)\g"(l2) >< 3(U)\3(Vi) > 



3(1 J) 



-1+1 



1 a n 



x (-1)^ I . JJ aj " 2 9J ~ 1 } TT V'a 

v 71 — J 

JV 

+ [ji en + 1) + j j(j j + 1)] n « s » on 

n=2 



(41) 



3.4.2 



7=1, J = 2 



E »(•") ( + 1 ) < 9(/J)l9"(12) X 9(/J)|a(12) > 



JV 

= a 2 (a 2 + 1) Yl S 9n a n 

n=2 

(42) 
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3.4.3 I = 2, J arbitrary 



E 92(1 J) ( g 2 (IJ) + 1 ) < g(IJ)\g" (12) X ff(/J)|a(12) > 



E 92(/J) ( 92 (/J) + l) 



(-i)«-^ +a>9S r +a )|.^ j£ J£ 



(-!)-<" A( ff/ ,aj) 



jj 97 
ji-i 97+ 1 ai 



(-^A( ai+1 , ai+l) {^ ; /+2 -+ + ;} 

(-^A( 9l+2 , ai+2) {^ 



X {-iyj-iA(g J _ 1 ,g^_ 1 ) 



3.1 9 j— 2 9.J 

jj-i gj g" 



"-iJ 



(-l)-°J-iA( Sj _ 1 ,a J _i) 



a.7-2 9.7 
jj-i 9.7 a.7 



::} 



I(-l)^-^-^+l-E^^n X(j2ijj) i A(a2) ^)|jl £ gj JJ A( S ;',a„)(-l)<- 1 +'»n- 1 + 1 {^ M 

^ ^ n— 3 ^ ^ 

f • 1 N 

V 7 \ JJ-l 9.7 a .7-1 J J- 1 , 9 " " 

n— J 

AT 

+ [j 2 (j 2 + 1) + jj(j.j + 1)] n 



3.4.4 



1 = 2, J = 3 



E 92(1 J) ( 92(/ J) + 1 ) < 9(/J)|r (12) X 9("M12) > 

= E 92(H) ( ff2 (/j) + 1 ) x 



x (• 



K-l^-M-iixO-a.jjji^aa,^')!^ 1 j f } (~^ { ? j f } 

JV 

+ [j 2 (j 2 + 1) + j 3 (j3 + 1)] I! 5 9>n 



(43) 



(44) 
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3.5 Explicit Formula for the Matrix Elements 

After having finished the precalculations in the last section we are now in the position to evaluate the whole 
matrix element in 1141 by using 1401 . 

We have now (again we abbreviate A(x,y) — \J (2x + 1) (2y + 1)): 
< a(l2)\q IJK \a'(l2) >= 

= E E 92(IJ)(g2(U) + l) < g(U)\g" (12) >< g(U)\S(12) > x 

ff" (12)9(7 J) 

x J2 92(JK)(g 2 (JK) + l) < g(JK)\g"(12) >< g(JK)\a(12) > 
g(JK) 

- [a(12) a'(12)] 



H3 



E 

g" (12) 



1 . ^i-i . ^j-i 



<*n, (_1)«J-1 (_l)-E„ = iJ™(_l) + l x 

x X(i / ,i J )iA(^',o / ){ a/ - 1 7 J} n A(s", „)(-1) 

, . „ v J-l N -| r JV 

x ( - i)aj { Y aZi 9j f/ \ n v„'-n n + n 



I 1 In dn— 1 



n — 2 n — J 



1 (_1JJK-JJ-E»=J I»-E* = itl Jm (_l) a .7-l (_1)»J-Sj (— 1)~ En=l 3« (— 1) H 



X {-ifK 

1 <%-l JK 



[g(12) <i=> a' (12)] 
2(-!) 



XO'^i^lA^J,^) a '-i a / [i J] A( fl '',a;)(-l)«n- I +<-x+ 1 P» 9 »ri f« 

J-l JV JV 

n 5 « n s « + [jjO'j+ i )+jifO'AT+ 1 )] n 5 « 



n=2 n=K 



if— lV J ^ ^ n=l J ™~Em = I+l •?>" (_l) a J-l (_l)aj-aj (— 1)~ En=l J" (_!)+! ; 



x X(j I ,jj)iA(a' I ,a I )\ ai - 1 £ 7J A«,a„)(-1) 



n=J+l 



^"n-l+ a »»-l + 1 J a « — 1 a ™ 

1 On On— 1 



x (-I)-" 7 JJ aj - 1 

1 1 aj-l JJ 



l( — iy K ~ JJ ~^n = l Jn-Z^ m = J+ i In. ( _ a j _ 1 ( _ a , - aj (_1)-E ra= i in (-1) + 1 X 



n=J + l 
J-l JV 

n=2 n=Jf 



-(— l) 3 '" JI "Sn=l Jn-E m=J+ i (_l)aj — Oj(-_]^- E„=i Jrt(_ 1)H 



XOV^UK.o/)] "Y 1 ^ ^ | n A«,a„)(-1) 



j-i 



(_1)«J J aj a J-l 
11 a J-l jj 



n=I+l 

jj(?j + i)+ijf(j'K + i) 



- | a^_ 1 +a„_i+l J J™ a n-l a n 

1 a n a n _i 



1-1 JV 

n—2 n—J 
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(-l) jK - j J-^'n = l i»-E«=J+l trn (_ 1 )a J _ 1 ^_l^a[,-aj (_^-E„=l Jn (-1) + ! ; 
X (-l) a K 



a K Jk Ofc-1 
1 a ' K -i 3 k 



n=J+l 

3l(jl + l)+3j(jj + 1) 



J-l AT 

II "W^, II S an.a'„ 
n=2 n = K 



+ i)+ix(iK + 1) 



n 



5(12) <^=> a' (12) 

(45) 

Here we have in the last step carried out the summation over g'{ 2 by evaluating all the ^-expressions. Finally 
we take a closer look at the symmetry properties which the four terms of the sum in I45H have with respect to 
the interchange [o(12) <=> a'(12)], that is the simultaneous replacement a n — > a' n , a n — > a n for all n = 1 ... N. 

• Fourth term : Due to the product of 5 ana > -expressions this term is obviously symmetric under 
[o(12) a (12)] 

• Third term : The symmetry is not obvious, we will show it part by part: 

1) In the (— l)-exponents aj-i = aj_i, a' K = cik by the (^-expressions at the end of the term. 

2) The term (— l) a J~ aj is not changed by interchanging a'j <-> aj, since a'j — aj is an integer number and 
therefore the first formula of 112111 holds. The integer-statement is verified by the fact, that S12, a' 12 are 
Standard Recoupling Schemes, defined as in l|A.2fl . Therefore they recouple all ji , . . . , jjv successively 
together according to theorem (IA. H . Since 012,012 contain temporarily recoupled angular momenta, 
namely a k = a k (a k - 1 ,j K ),a' k = a' k (a' k _ 1 , j k ) for which \a k -i - jk\ < a k < a fc _i + j k ,\a' k _ 1 - j k \ < 
a ' k ^ a 'fc-i + jk the integer- or half- integer- property of each component a k ,a' K is only caused by 
the order the involved spins ji, . . ■ ,Jn are coupled together. Since this order is the same in 012, a" 12 , 
the components a k ,a' K are simultaneously (for every k — 1 . . . N) either half- integer or integer and 
therefore every sum or difference a k ± a' k is integer. 

3) The same statement as in 1) holds for a' J _ 1 ,a' K as the entries in the upper left corner of the two 
6j-symbols before and after the product in the middle of them. 

4) In the product of the 6j-symbols the exponent of the (—1) contains only a sum a n _i + a' n _ 1 and is 
therefore symmetric. 

5) All prefactors A(a k , a' k ), k = J . . . K — 1 are symmetric, too. 

6) Finally all 6j-symbols in the third term turn out beeing symmetric, if we recall, that they are invariant 
under an interchange of their last two columns l|122|l followed by an interchange of the upper and 
lower arguments of their last two columns. 11231 . 

• Second term : The symmetry is again not obvious, so we will show it part by part: 

1) In the (— l)-exponents a/-i = a' I _ 1 , a'j — aj by the 5-expressions at the end of the term. 

2) The term (— l) ai ~ a i is again not changed by interchanging ar <-» a'i, by the same integer statement 
as under point 2) in the Third term discussion. 

3) The same statement as in 1) holds for a/_i,aj as the entries in the upper left corner of the two 
6j-symbols before and after the product in the middle of them. 

4) In the product of the 6j-symbols the exponent of the (—1) contains only a sum a' n _ 1 + a n _i and is 
therefore symmetric. 

5) All prefactors A(a' k , a k ), k = I . . . J — 1 are symmetric, too. 

6) Finally (again) all 6j-symbols in the third term turn out being symmetric, if we recall, that they are 
invariant under an interchange of their last two columns l|122|l followed by an interchange of the upper 
and lower arguments of their last two columns. 11231 . 
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• First term: This term is not symmetric under [a(12) a' (12)]. Arguments similar to the ones we 
gave in the previous points now let us conclude that: 

1) The prelactors (— l) aj , (— l)"- 7 - 1 are not symmetric with respect to the interchange ol 
a.j — ► a'j,a' J _ 1 — > aj-i. 

2) The two 6ji-symbols containing a,j,a' J _ 1 in the upper left corner are not symmetric either. 

The analysis has revealed that the last three terms in the sum of {HJ are symmetric in a(12) and a' (12), 
hence after an antisymmetrization with respect to the interchange of a(12) and a' (12) only the first term in 145H 
survives. 

Summarizing, we get by explicitly writing down all the terms occurring through the antisymmetrization: 
< a(12)\q IJK \a'(12) > = 



1 a n a 



n Q>n— 1 



1 1 aj-l 3J 



(_1)JA'-J./-E„ = i J»-E m= J + l J-"(_l) a J-l(_l)a' I -a I EjUl in(_l) + l x 



x X(Jj,j K )2 A(a'j,aj) 



1 jj 



X (-l) a K 



1 <_i hi 



J J n A(a„,a^)(-l) a »-^ 

J ' n=J+l 

1-1 AT 

n=2 n = 



l(_l)J-'-^-E^ = i 3«-E TO= / + i ?>"(_l)<4-l(_l)aj-ar(_l)-El=l J"(_1) + 1 X 

I 1 J./ 



(— l) 3K_JJ_ E ra= i Jn-E m=J + i Jm (_j)a' f -af £„=J in (— l)" 1 



! 1 ajc_i ok 



n = J + l 
1-1 AT 

n=2 n = 



1 J Jn a„ 
1 a n a n -i 



(46) 

Here we have underlined the terms, which are different with respect to the antisymmetrization (recall, that 
(— l) ai ~ a i = (_i)-( a /- a /) = because the exponent is an integer number. All the othe r term s are 

symmetric under the interchange a(12) a'(12), again because of the symmetry properties 1)1220 . I|12M0 and 
the symmetrization by the 5-expressions. 

Before we write down the final result we can simplify the exponents in l|46[l : 

EI — 1 ^J— 1 . v-, J 1 ■ *^J — 1- v-^/f — 1 „7_i 

( _j , n = lJ«-l-m = I+l^(_l)-E Tl= iJr l (_ 1 )-l-n = l^-^m = J+l^(_l)-E ri = i 3n = 



= (-i)- 2 E„ = iJ n (_i)-E„ =J+1 j n( _ 1) -2E„ =1 (_ 1 )-E„ =J+1 ^ 
= ^^=/+i : '"(_i)- 2 E; l= j^(_ 1 )- ^„=j+i ^ 

= (-l) 2 ^+E„=/ + l J"-E„ =J+ i 3n 



23 



Now we are able to give a closed expression of the matrix elements of qijK in terms of standard recoupling 
scheme basis 1141 . In order to avoid confusion we assume that / > 1, J > I + 1, the remaining cases will be 
discussed below: 



THEOREM 



< a\qij K \a > = 



xX(j I ,jj)^X(jj,j K )^(2a I + 1)(2< + l)yj{2aj + l)(2a' J + 1) x 



ai-i Ji ai 
1 Oj ji 



JJ V(2< + l)(2a„ + l)(-l)<-i+^-i+ 1 J J " a " 

n=/+l v 



TT V(2< + l)(2a„ + *• a ™- 1 a " 

n=J+l v 



1 1 aj-i JJ J 1 1 aj jj 



(-1) 



1 a' A '-i 3K 

aj+aj-l J a J a J-l 1 J J J a J 



1 a.j-i jj 



1 a.j jj 



x JJ <5 a „ a ; JJ 5a <>< 



(47) 



with X{ji,j2) = 2j'i(2ji + l)(2ji + 2)2j2(2j2 + l)(2j2 + 2). Notice that all still appearing 6j-symbols are just 
abbreviations for the following simple expressions in which no summations or producs (factorials) need to be 
carried out any longer as compared to (II 20t e.g. (using s — a + b + c) : 



a b c 
1 c b 



a b c 
1 c- 1 b 



a b c 
1 c- 1 6-1 

a b c 
1 c- 1 6+1 



= (-ir +1 
= (-ir 

= c-ir 



2[6(6 + l)c(e+l)-a(a+l)] 
[26(26+ l)(26 + 2)2c(2c + l)(2c + 2)]3 

2(s + l)(s - 2a) (s - 26) (s - 2c + 1) 
26(26 + 1)(26 + 2) (2c - l)2c(2c + 1) 

s(s+ l)(s - 2a - l)(s - 2a) 
(26 - 1)26(26 + l)(2c - l)2c(2c + 1) 

(s - 26 - l)(s - 26) (s - 2c + l)(s - 2c + 2) 
(26 + 1)(26 + 2) (2b + 3) (2c - l)2c(2c + 1) 



(48) 
(49) 

(50) 

(51) 



Remark: Gauge invariance 

Recall that a vertex is said to be gauge invariant if the angular momenta coming from the edges ei, . . . , ejv 
meeting in the vertex v are coupled to a resulting angular momentum j = qn — 0. Using the notation from 
1A.5II this means: 



Gn-i + Jjv = 



This implies: 



Gn-1 = GjV-2 + Jn-i = — Jn 
~> GjV-2 = GjV-3 + JjV-2 = — JiV — Jn-1 



(52) 



(53) 
(54) 



But that gives gjv-i = Jn and a certain restriction on which values <?at-2 can take due to the Clebsh-Gordon- 
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theorem <|A. ip : 

max(|jjv_ 2 - ffjv-sl, b'jv - JJV— 1|) < 9N-2 < min (jjv-2 + gN-3 , Jjv + jiv-i) (55) 
This relation will become useful when we consider gauge invariance later. 
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As promised we now display the remaining cases of l|47^ explicitly. They are obtained, if some of the special 
cases I4H . . . 14411 are involved. 



3.5.1 1=1 J,K arbitrary 



<a\q 1JK \a > = i(-l)^-«+°^+ 1 (-l)E»=2^(-l)-Sf = j 1 + iJ !1 x 



xX(j u jj)*X(jj,j K )*^paj + l)(2o'j + 1) x 

J-i f • ' ' 

jj v(2< + i)(2a„ + i)(-i)"»-i+-»-i+ i J ; a "- 1 °» 

„ I 1 tin CSn-1 

n = 2 v 



— l+ a n — 1+1 J Jn a n — l a n 

1 a n a n _i 



JJ V(2< + l)(2a n + l)(-l) 

n=J+l 



1 a'if-i 



f_ifj+«j-i J aj J J a ^-i I J a J-i J J a J I _ J a J a J-i [ J a J-i 3-1 a J 

1 1 l 1 aj_! jj J 1 1 aj jj J 1 ; 1 aj-x jj J 1 1 aj jj 



3.5.2 1=1 J=2 K arbitrary 



<a|<?i2A-|« > = i(-l) JA "-^--' 2+ ' lK+1 (-l) E ;.=2^(_l)-S P= .7 + iJp 



xX( j2 ,j if )V(2a 2 + l)(2^ + l){ J ] 1 g ^ 2 } 
JJ V(2< + l)(2a n + l)(-l) 

A? 

n* 



^a^_ 1 +o Tl _i+l J a n— 1 tt n 

1 a n a n -i 



ClK JK a-K-i 
1 ljf-1 JK 



02(02 + 1) - a 2 (o2 + 1) 



(57) 



3.5.3 



1=1 J=2 K=3 



This case is actually the easiest. We will use it in the next section. Therefore we will write down the calculation 
explicitly. We start with JI3J to obtain 



< a|<?i23|a' > = a 2 (a 2 + 1) <5<'a„ 



l(-l) h - jl+j3 X(j 2 ,js)i^(M + l)(2a> 2 + 1) { j l * f\x 

x (-ir a l J J IJ + + !) +m*> + 1)] IJ <w< 



a < > a 



= a 2 (a2 + 1) - a 2 («2 + 1) 



^(-lY 2 - J1+33 X(j2,h) V(24 + l)(2aa + 1) | * J? J | x 

AT 

IJ *»»« ; ( 58 ) 



x (-1)" 3 ^ J * " J 



«2 J3 



(56) 



2G 



Now the conditions for the arguments in the definition of the 6j-symbols l|12()fl gives certain relations for a 2 , 
a' 2 occurring in 1581 . namely: 

(0,-1 

a 2 - < a,2 

I 02 + 1 



a 2 = a,2 



We can choose either the first or the third case to obtain a non- vanishing matrix element. We choose 
(the other choice would give us a sign, due to the antisymmetry of q). So we continue (considering only the 
nontrivial information a 2 , a' 2 contained in the recoupling-schemes a, a' ): 



< a 2 \q 123 \a 2 - 1 > = a 2 (a 2 + 1) - (a 2 - l)a 2 



±(-l)'»-*+*XCfr,is) V(2(oa - 1) + l)(2oa + 1) { * ja _ ± " 2 



x (-ir ? J:i i a . 2 

{ 1 ai - 1 j 3 

Now we can rewrite the second 6j-symbol in l|59|l by using the identity 2 



a b c 
1 c- 1 b 



2a(2a + l)(2a + 2) / b a c 



2b(2b+ 1)(26+ 2) [lc-la 



< a2|gi23 a2 - 1 > 



2a 2 ^(-l)^ 1+j3 X(j 2 , j 3 ) V(2a 2 - l)(2a 2 + 1) 



Jl J2 



«2 



1 a 2 — 1 j 2 



x (-!)»» 



2a 3 (2a3 + l)(2a3 + 2) 
2j3(2j 3 + l)(2j 3 +2) 



J3 «3 £12 
1 a 2 — 1 03 



At the last step we once more use B5I to express all 6j-symbols in 1601 explicitly. Furthermore we use the 
fact, that ]3 + a 2 + 03 is an integer number and therefore (_i) 2 C?3+ a 2+a3) = 1 (one can see this by applying 
the integer conditions 1B.2II to the second 6j-symbol in I60H ). Additionally remember the shortcut introduced 
earlier: X(j 2 ,jz) = 2j 2 (2j 2 + l)(2j 2 + 2)2j3(2js + l)(2j3 + 2). After carefully expanding all the terms and 
cancelling all identical terms in the nominator and the denominator the explicit result is: 



(t 



< a 2 |<2i 23 |a 2 - 1 > = 



V'(2a 2 -l)(2a 2 + l) 



(,7i + h +a,2 + 1) (- ji + h + a 2 )(ji - ji + a 2 )(ji + ji - a 2 + 1) 



(js + a 3 + a 2 + 1)(-J3 + a 3 + a 2 )(j 3 - a 3 + a 2 )(j 3 + a 3 - a 2 + 1) 



= - < a 2 — l|gi 23 |a 2 > 

(61) 

The analytical result for this special case coincides with the result already obtained by graphical methods 
in |S], if one considers the gauge invariant 4-vertex (see next section!) by putting a 2 — > ji 2 and a 3 — » j\ (as a 
consequence of applying the definition of the standard recoupling schemes (definition IA.20 to a 4-vertex). 



3.5.4 Comparison: Computational Effort 

At the end of this section we want to compare the computational effort one has to invest for calculating the 
matrix element 1141 using the full definition in terms of 6j-symbols I20H or the derived formula 14711 instead. 
We will give here only a rough estimate, since for the full definition 1201 the calculation can hardly be done for 
all possible combinations of arguments. 

Consider an A-valent monochromatic vertex 3 v with iV outgoing edges e%, . . . , ejv, each carrying the spin 
ji = . . . — j'jv = jmax- Assume, we had to calculate the matrix element (|14^ for a certain combination of the 
triple I < J <K, namely / « (J - 1) « (K - J) := L ~ f > 1. 

2 This follows directly from 1491 

3 Of course this special case is not that expensive, because it is the most symmetric case. But it illustrates the estimates for the 
general case with different spins. 
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Let us first discuss the full definition using il I It with (I2( H inserted. 
Consider first the definition 11201 of the 6j-symbols: we will only pay attention to the ^-coefficient, since the 
number of A-coefficients is constant. By the requirement for the summation variable n (max[ji + 32 + j'12, ji + 
j + 323 , 33 + 32 + j'23 , 33 + j + J12] < n < min[ji + j 2 + 33 + j, h + J12 + j + hi , J12 + ji + 323 + js\) we can (if 
we put ji — 32 — js = j = jmax, 323 = ji2 = 2j) extract < n < 4j. That is we have approximately Aj max 
summations and therefore 7 • 4jmax factorials to calculate for every 6j-symbol. 

Now look at the definition of the 3nj-symbol in terms of 6j-symbols 1181 : We have approximately 21 — 2L 
6j-symbols, due to summation over the intermediate recoupling steps hk and additionally J — I = L 6j-symbols 
not involved into that summation, a constant number which we can drop. Now in the worst case: 

h\ ^ 2jmax , 3 max 

< ti2 < 2j max + jmax, - ••,(-?- l)jmax < hi-i < (I + l)j ma x- Therefore each h k 
(2 < k < (I — 1)) can take 2j ma x different values and we have thus about (2j m ax) 1 possible combinations for 
the hk. So we had to calculate at all 2I(2j ma x) 1 = 2L(2j ma x) L 6j-symbols. 

Every term in the sum (lilt contains a product of 4 3nj'-symbols. Now the summation over g(IJ), g(JK), g" (12) 
again gives (under the assumption, that each intermediate angular momentum g(IJ)k, g(JK)k, g(12)' k ' can take 
2jmax different values, but only ( J — I) « (K — J) := L intermediate steps of each summation contribute to 
calculations, due to the 5-terms in I I fifi ), for each matrix element approximately (2 • 2j ma x) iL 3nj-symbols to 
calculate. 



Summarizing these 3 steps we end up with a computational effort of approximately: 

7 " 4j maa; ■ 2L(2j ma x) ' (2 ' 2j rna x*) ^ (jmax) ^ (jmax) ^ 

calculations of factorials occurring in l|120|l . 

It is much easier to discuss the effort one has in case of using the derived equation l|47p. We have only a 
product of special 6j-symbols containing no summation and factorials at all. So one only has to carry out 

the product consisting of only 



K-I = 2L~ fN factors 



, independent of j r , 



Be aware that this estimate given is only rough, one could introduce the symmetry properties l|122|l . I|123|l 
and additionally keep in memory previously calculated 6j- or 3n7'-symbols to save calculation time. Nevertheless 
the computational effort for the calculation of the matrix element 1)14^ depends on jmax if one uses the original 
formulas JT]] with PD). This is no longer the case if one uses 147H . It is clear, that if one wants to numerically 
compute all the matrix elements then one cannot get very much over jmax = 2 with 1141 . 



3.5.5 Conclusion 

We have shown in the last section, that it is possible to explicitly evaluate the matrix elements of the volume 
operator in l|14fl . Here by 'explicitly' we mean that there are no more 6j'-symbols in the final expression. 
The derived formula is a simple algebraic function of the spin quantum numbers, no factorials appear any 
longer and no conditional summations, implicit in Racah's formula for the 6j-symbol, have to be carried out 
anymore. Thus the computational effort in order to evaluate the matrix elements has decreased by a huge 
order of magnitude, which grows with growing maximal spin jmax- This simplification has been achieved by 
the discovery of a nontrivial fact, namely that the highly involved formula 114H or 120H is like a telescopic sum 
of the form X^=i( a « — <*n-i) = a Jv — ao once one takes the orthogonality relations of the 6j-symbols and the 
Elliot-Biedenharn identity into account. 

A first observation is that the matrices defined by l|47|l . show a banded structure, that is a rich selection rule 
structure. Non-vanishing entries are only on certain parallels to the main diagonal, because of the restrictions 
of the presence of an entry 1 in every 6j-symbol contained in 1471 . 
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4 Gauge Invariant 4- Vertex 

In this section we will examine in more detail the gauge invariant 4-vertex, that is the following configuration 
of edges: 




Figure 1: The configuration at the 4-vertex: 4 out- 
going edges each carrying a representation of SU(2) 
with a weight according to ji, j 2 , h, n 



We have 4 edges ei, . . . , e$ outgoing at the vertex v carrying the spins ji, . . ■ , ji and the according represen- 
tations Of SU (2), TTj 1 , . . . , 7T J4 . 

The square Q v = V% of the volume operator represented in term of the standard-recoupling scheme basis 
was (note that in what follows we will stick to the squared version of the volume operator, therefore of all 
eigenvalues we write down the square root has to be taken in order to obtain the spectral behaviour of the 
volume operator itself): 

Qv ■- z ■ < L .J> K ) [(Ju) 2 ,(Jjk) 2 ] 

KJ<K 

= Z- Y, < L . J > K ) fax ( 62 ) 

KJ<K 

Since we have 4 edges the summation in 162H has to be extended over the combinations: 
(KJ<K) = (1, 2, 3), (1, 2, 4), (1, 3, 4), (2, 3, 4). 

The point is now that due to gauge invariance the 4 angular momenta j\ , . . . , ji should couple to a resulting 
zero angular momentum j = at the vertex v. For this reason for the angular momentum operators Ji, . . . , J a 
holds due to l|55|l : 

Ji + J 2 + Js + Ji = (63) 

which implies 

Ji = -(Ji + J 2 + Js) (64) 

It follows that: 



9124 = — [Jll'Zl + 9122 + <?123 j = —9123 
9134 = -f9l31 + 9132 + 9133) = -9132 

9234 = — (9231 + 9232 + g233j = —9231 (65) 



Here we have used the fact that qijj + qui = V/, J 4 . 

4 To see this, just take the definition qijK = |^(^/j) 2 i (Jjk ) 2 j an d expand the resulting commutators. Alternatively use the 
antisymmetry of qijj ~ tijkJjJjJ 1 } to see that qijj + qui = quj — qjn = . 
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Thus l|62[l reduces to: 

e(I, J, K)q IJK = e(l, 2, 3)<?i 23 + e(l, 2, 4)q 124 + e(l, 3, 4)<?i 34 + e(2, 3, 4)q 234 

KJ<K 

= [e(l, 2, 3) - c(l, 2, 4) + c(l, 3, 4) - e(2, 3, 4)] g 123 

= 2 ■ $123 



(66) 



where we have used in the last line the configuration of the 4 edges outgoing from v described above to obtain 
the orientation (±1) of every triple of tangent vectors corresponding to 3 distinct edges (Note that for different 
orientations of ei, . . . , 64 just the prefactor 2 changes.). This brings us precisely into the situation of the last 
example in the previous section. We can now explicitly write down the matrix-elements of Q v represented in a 
basis of standard gauge invariant recoupling schemes |a(12) j = M — >= |a 2 >, \a'(12) j = M = >= 
|a 2 >. For simplicity we relabel a 2 — + j\ 2 , a 3 = a' s — *■ 34, a 2 — » j' 12 , o 3 = 0,3 - * J4- 
Now the non vanishing matrix-elements in 1)61^ are: 

< ji 2 |<7i 23 |ji 2 - 1 > = 



y/(2fa - 1)(2 J12 + 1) 



O'l + J2 + jl2 + l)(-jl + J2 + jl 2 )(jl - J2 + jl 2 )(jl + J2 - jl2 + 1) 



(j 3 + J4 + Jl2 + 1)(-J3 + J4 + jl 2 )(j 3 - J4 + jl 2 )(j 3 + J4 ~ jl2 + 1) 



< jl 2 - l|<?l 23 jjl 2 > 



(67) 
(68) 



By 15511 we get certain restrictions for the values that ji 2 may take 5 : 

max (| ji - j 2 \, \js - 3a\) < ji2 < min (ji + 32,33 + ji) 
Therefore the dimension n of the matrix-representation A of <Ji 23 in the standard basis is given by. 

n := dim A = min (ji + j 2 , 33 + ji) - max - j 2 \, \ja - h\) + 1 

• max -min . -1 frr\\ 

= 3i2 -J12 +1 (69) 

We find for the matrix A (labelling the rows by ji 2 and the columns by j[ 2 , where the first row/column equals 
ji2 = j'12 = 3r2 n increasing down to the last row/column with j i2 = j' 12 = jvf x an d using the abbreviation for 
the matrix element au '■= i' < jii" 1 + fc|<?i 23 | ji2 n + k — 1 > (where i is the imaginary unit 6 ), k = 1, ... ,n — 1 
and ao = a n = 7 ): 



A = 



( -ai 
ai — a 2 
a 2 



V 



\ 



— a n -i 
a n -a / 



(70) 



That is: The matrix A possesses a banded matrix structure which is called a Jacofei-matrix. Note, that the 
0^ are purely imaginary, because A is hermitian and its eigenvalues are real. We will discuss the spectral theory 
of A by analytical and numerical methods. Note the following advantages of the gauge invariant case over the 
gauge variant: 

• The dimension of A scales only linearly with the spins outgoing at the vertex v (this advantage will be 
useful for the numerical studies). 

• There is no sum over matrices left any longer, as it would be the case for the gauge variant 4-vertex. 

• Due to the formulation in a recoupling scheme basis we have automatically implemented gauge invariance. 



J We may label w.l.g. edges in such a way that < j\ < j 2 < j 3 < j'4 < (ji + j 2 + J 3 ) and ji + j 2 + jr 3 + j'4 =integral 
6 This only changes the spectrum of A from being antisymmetric to hermitian and therefore rotates its spectrum from purely 
imaginary to purely real. 

7 just insert j\ 2 = J™" or J 12 = ill" 1 + n = 3i 2 ax into the matrix element 1671 
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4.1 Analytical Investigations 
4.1.1 Eigenvalues 

As pointed out before, all eigenvalues A of A are real and come in pairs ±A. The special case of zero-eigenvalues 
will be discussed below. One can find upper bounds for the eigenvalues by applying the theorem of Gersgorin 
(see |Hl,p.465): 

Theorem 4.1 Gersgorin 

Every characteristic root A of a (n x n)-matrix A lies at least in one of the discs 



\au — A| < 



i = 1, 



That is, every eigenvalue lies in a disc centered at the diagonal element an with radius of the sum of moduli 
of the off-diagonal-elements aij, i ^ j of the i-th row or column (called the i-th row- or column-sum). In case 
of the gauge invariant 4-vertex this theorem simplifies due to the banded matrix structure of 17UH and the fact 
that an — to 

\M < ^2 = \a,i + ai i+ i (71) 

We will give an upper and a lower bound for the eigenvalue-spectrum in terms of the leading polynomial 
order of the largest angular momentum j ma x ~ max(ji, . . . , j^). 

By inspection of 1)70^ one can see that the row sum introduced in theorem i|4. 1 pi is dependent of the modulus of 
each of the matrix elements ai(ji 2 ) :=< ji2 | qwi \ ji2 — 1 >, i = l...n— 1. This observation will be useful 
for obtaining an upper bound for the modulus of the eigenvalues of A. On the other hand we could also use 
theorem l|4.1[l for giving a lower bound of the eigenvalues if we could guarantee the existence of the inverse A -1 
that is the absence of the eigenvalue 0, which will be discussed explicitly in appendix[U] Then the upper bound 
of the eigenvalues of A" 1 would give us a lower bound of the non-zero eigenvalues of A. However, due to the 
general formula for the matrix element of the inverse matrix 



(A hj -^A~ 



(72) 



(where Mtj denotes the sub-determinant of A with row i and column j deleted) all the entries of A 1 will 
be of the order a .(j 12 \ ■ We can therefore try to find the extrema of the matrixelement Hti7fl in terms of the spins 

ji, . . . ,j4 by partial differentation 8 . Note, that we have the freedom to choose ji < ji < 33 < j'4 = jmax, since 
the matrix element 167H is symmetric under permutations of ji, . . . , j'4 : 



d 


a(jl2) 




dji 


d 


a(ji2) 




dj 2 


() 


aO'12) 




dji 


d 


aO'12) 



O j : 



(o) 



djt 



(o) 



= ji 



1 1 

2 + 2 
1 1 

-2 + 2 

1 1 

2 + 2 

1 1 

2 + 2 



1+4^2 


+ 4.72 


+ 4jI 


l + 4j? 2 


+ 4ji 


+ 4j? 


l + 4j? 2 


+ 4J4 


+ 4j 4 2 


l + 4jr 2 2 


+ 4js 


+ 4j 3 2 



(73) 



If we want all relations of 1731 to be fulfilled at the same time with strictly positive values of ji , . . . , j'4 then 
we have to demand 



ji = 32 ■■= m and j 3 = j 4 := jmax 
Therefore l|69^ reads due to the ordering j\ < j 2 < 33 < J4 = jmax as 

< ji2 < 2m 

And the matrix element l|67|l simplifies to 



(74) 



(75) 



a(ji 



i 2 2 



v 7 ^ 



2 - 1 
12 L 



[(2a+l) 2 -j 2 2 ] [(2jmax + i f ~ j 



(76) 



3 We will list here only these solutions which allow positive values for Ji, . . . ,J4 and mutually different ji, 32, 33,34, 
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1. Largest Eigenvalue 

By inspection of 1761 we can maximize the order of jmax contained in the matrix element by putting 

m ~ jmax ~> a(j 12 ) ~ jm ax (77) 

The result is an upper bound on the growth of the maximum eigenvalues of the matrix A with the maximum 
angular momentum jmax 9 : 

3 3 

| Arnax (jmax ) | ^ jmax = ^ I Vmax (jmax ) | ~ jmax (^^) 



2. Smallest non-zero Eigenvalue 

By assuming the existence of the inverse A -1 with a sparse population of entries of the order — one can 
minimize the order of jmax contained in the matrix element (that is maximizing the matrix elements of 
j4 _1 ) by putting 

171 ~ 1 <S jmax ~> a(j 12 ) ~ jmaz (79) 

The result is an upper bound on the growth of the maximum eigenvalues of the matrix A -1 with the 
maximum angular momentum jmax and therefore for the smallest non-zero eigenvalue of A 10 : 

i 

| ^min (jmax ) | ~ jmax | Vmin (jmax ) | ~ jmax V^^V 



These are first estimates, we will come back to this, when we discuss the numerical investigations, since the 
theorem (14. II and our approximations tell us nothing about the numerical coefficints in front of the leading 
powers of jmax- Nevertheless this estimate will give us a certain criterion for completeness of numerically 
calculated eigenvalues: Since the smallest eigenvalue \ m in grows with j mM we can at a certain value of jmax be 
sure having calculated the complete volume-spectrum for all V < V m in as we shall see in the numerical section 
below. 



4.1.2 Eigenvectors for A = 

Posing the eigenvalue problem = A4< for the matrix A we obtain a three term recursion relation every 
eigenvector $ of A has to fulfill: 

at-i^k-i — fflfc^fe+i = A^fc with ao = a n = 0, A 6 M, ak(0 < k < n) purely imaginary (81) 

We can now check, whether the eigenvalue A = belongs to the spectrum. This decouples the recursion 
relation 1811 to give 

ak-i^k-i ~ «fc*fc+i = (82) 

Now for consistency of l|82l : 



k=n: 



k=l: 



o n -i*n-i = and therefore * n _i = ^ n - 3 = . . . = 

-cii*2=0 and therefore * 2 = *4 =...=0 



But this means that the matrix A can only have an eigenvector ^= if the dimension n of A is odd, because 
if n would be even, then all components of $ would be forced to vanish by l|83^ . 

That means, that we will only obtain A = as a eigenvalue in configurations with odd dimension of A. We 
can now construct explicitly the eigenvector ^ for odd n: First we choose <&i = x with x — const 6 C Then 
from 

O.k-1 ^fc + l 



We find the general expression: 



(84) 



n a 2s-i 

*2r+i = x ■ where r = 1, 2, . . . , a=i (85) 

Yl 0,2s n = dim A 

s=1 x = *1 



9 This result coincides with that already obtained in [241 . Note, that the number of terms in each row/column - sum is equal to 2 
due to the special structure 1701 of the matrix A and therefore independent of j m ax ■ 

10 Note, that one should really use (A~ 1 )ij which can be a complicated polynomial in the Aij. This is indeed the case as shown in 
appcndixE]and therefore the result presented here is at best a rough estimate. 
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Since the a's were chosen to be purely imaginary, all odd components $3, . . . , 'fn of 9 are real, all even 
components ^2, • ■ • , ^n-i are identical zero, because of 1831 . Finally we can fix ^1 = x to be (up to the sign of 
x = #1): 



X = *1 = ± 





r 


2-1 


n-1 

2 


n a 2 S -i 




i + E 


s = l 




r 






n 

s = l 





(86) 



Summarizing, the eigenvalue A = only occurs in the spectra of matrices A, possessing an odd dimension 
n, as a single eigenvalue (since its eigenspace is only one-dimensional due to the uniqueness of construction of 
the eigenvector $ up to a constant rescaling). Hence we have shown that A is singular iff n odd and in that 
case A = has multiplicity 1. 

4.1.3 Monochromatic 4- Vertex (ji = fa = j'3 = j'4 = j) 

Observing this special case, the matrix elements in (jri7^ simplify dramatically to 



< J12 |<7l23 I J12 - 1 > 
where < j'12 < 2j and dim A — n — 2j + 1 



- 1 



0i 



(87) 
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4.2 Numerical Investigations 

In this section we will describe numerical calculations done for the gauge invariant 4-vertex. We will (after 
describing the setup) sketch the computational effort first. Secondly we will give a conjecture about the volume 
gap, that is the smallest non- vanishing eigenvalue as a result of the calculations. As a third step we will take a 
look on the accuracy of the upper bound given by theorem l4.ll Finally we will present some spectral estimates. 



4.2.1 General Setup 

We calculated for the gauge invariant 4-vertex the spectra of all possible edge-spin-configurations ji, j 2 , j3, ji 
up to a maximal spin of j m „ = 50 using the mathematical software Maple 7. Since the matrix element 167H is 
symmetric with respect to interchange of the j's we calculated the spectra of all Qv{ji,32,32,jz) for 

< ji < 32 < js < ji < min(ji + j 2 + ja,jmax) and j 1 + j 2 + 33 + 34, integral (88) 

Thus we compute less than ^ configurations, without losing any information. The conditions on the right side 
of 1881 ensure, that we exclude all trivial configurations, since if ji +32 +ja < ji or j\ + j 2 +js +34 not integral 
it would be impossible to recouple them to resulting zero-angular momentum, being the definition for gauge 
invariance. If one spin would be identical zero, then the obtained configuration would also be trivial, since 
it would describe then effectively a gauge invariant 3-vertex, which vanishes identically, hence we also impose 
3l ■ h ■ 33 ■ ji > 0. 

The possible values for the intermediate recoupled ji 2 are then due to l|68|l and the above introduced order 
of the angular momenta (since we have sorted the ji, . . ■ ,ji by their modulus and each 3 > 0, we can leave out 
the modulus-notation in I|ri8jl by writing the j's in certain order) : 



max(j 2 - j\,ji ~ 33) < 312 < min(ji + 32,33 + ji) = ji + .72 
The dimension of the matrix A of such a configuration is then according to 1691 given by: 

dim A = dim(ji , j 2 , 33 , ji) = j\ + 32 - max(j 2 - 31 , 34 -33) + ! 



(89) 



(90) 



For every configuration then the matrix elements according to 167H are calculated and inserted into a nu- 
merical matrix. This matrix is then numerically diagonalized, its eigenvalues are sorted ascending and from this 
spectrum all eigenvalues > are taken. These data are then written into a file linewise, each line starting with 
the values of j\ , . . . , ji and the total number of saved eigenvalues followed by the sorted list of the eigenvalues 
itself. Of course we have to keep in mind the multiplicity 2 of every saved eigenvalue > 0. Additionally we 
have to pay attention on the ordering procedure we applied on ji .. . ji whenever we work with the number of 
eigenvalues, since we have suppressed certain multiplicities. The following table gives the resulting multiplicity- 
factors by which eigenvalue numbers resulting from corresponding spin-configurations should be multiplied 1 : 





all spins different 


one pair equal 


two pairs equal 


three spins equal 


all spins equal 


no ordering 


TT =24 


#-12 


2!-2! ~ u 


3! * 


£ _ 1 

4! 1 



'In general we have N\ possibilities to arrange a list of elements, but having M identical elements each of them with multiplicity 



K 1 , K 2 , . . . , Km in that list we can only have 



JV! 

K 1 \-K 2 \-...-K M \ 



different arrangements of the elements of that list. 
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4.2.2 Computational Effort 

According to the setup above we would expect that: 



• Number of Configurations 

First of all we choose integer spins for simplicity, that is jj — > a.j — 2 • jj and assume all a's to be 
different (We therefore neglect configurations containig equal spins which are dominated by those with 
differnt spins). We choose one of the 4 integer-spins to be maximal and constant, say a a- = a max := 
max (at), L = 1 ... 4) and then by condition l|88[l we have clk =< min ( a L, o-max)- Therefore a max < 

Y aL- Labelling the three remaining aL by a\, 02,13 we get from that 03 > a max — 0,2 — 01. But this 

LjLK 

is only the question of counting the number N of points of a three dimensional cubic lattice fulfilling the 
last condition which is given by: 

a max a max a max a max — 2 drnax — a l a rnax — a 1 — a 2 

N(a max) =Y: E Ei- E E Ei (3D 

01 — 1 Q2 — 1 «3 = 1 ai=l (12 = 1 03=1 

Cmax 

Finally we add Yl N(a max ). The result is the total number of calculations ff(c m „). Due to the integer 

condition in 1881 . which becomes an even-number condition in the a's, we have to divide that number by 2 
and to multiply it by 4, since we have chosen one of the 4 integer spins to maximal arbitrarily. If we finally 
plug in Cmax = 2 ■ jmax, we get for the number of configurations N(j ma x > §) with all spins different: 

20 3^2 

(jmax^) = ' jmax "t - 12 - j max ~t~ ^ ' jmax 5 • jmax ~\~ 2 (92) 

This can be compared to the numerically fitted curve of the number of configurations N num '- 

N num (jmax) = 6.67 ■ jtax + 13.33 ■ j m ax + 11.42 • jmax - 8.30 ■ jmax + 2.22 (93) 




jmax 

Figure 2: Number of configurations N(j max ) 

• Number of Eigenvalues 

To calculate the expected number of eigenvalues E(j m ax we have to sum over the dimensions of the 
individual representation matrix of the Volume Operator which is given by 

dim = min (ji + j 2 , j 3 + ji) - max (| j 2 - ji | , \ja - 33 \) (94) 

One would expect that E(j max ) ~ jmax- From the numerical calculation we get for the total number of 
eigenvalues (including 0-eigenvalues and multiplicities): 

E(jmax) = 2.67 ■ jmax + 10.00 ■ j m ax + 15.52 ■ f max + 7.92 ■ j max + 27.90 • jmax - 77.23 (95) 
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jmax 

Figure 3: Number of eigenvalues obtained in 
dependence of j max 



4.2.3 First Impressions 

We drove the calculations up to a value of jmax = 50 12 . Our primary goal is to obtain some hints on a possibly 
analytical eigenvalue distribution function for large jmax- 

Therefore as a start we scanned through all the configurations saved in a file and computed an eigenvalue- 
density by distributing the possible eigenvalues A into discrete intervals of width AA = 0.5 where A belongs to 
the interval /„ = round(-^ + 0.5) = [A n — AA, A n ] where A n = n ■ AA, n = 1, 2, 

For this we take the square root of the calculated eigenvalues, since we computed the spectrum of the square 
of the eigenvalues of V v and the interval width AA does not scale linearly if we would first sort in the eigenvalues 
of Q v and then took the square root. Furthermore, we drop the prefactor Z in what follows. 

Then we plot the logarithm+1 of the total number of eigenvalues > (including their multiplicity) according 
to the multiplicity tabular given before) in a certain interval /„, against 2 • jmax and the volume-eigenvalues 
denoted by 2 • V (the form of the axes-labels with prefactors 2 as well as the +1 in the logarithm are only for 
technical reasons). The result is given in figure H"2.MI 

The diagramme suggests, that for each interval only configurations up to jmax matter. Configurations 
with jmax > jmlx do not increase the number of eigenvalues in the interval Therefore we are led to the 
idea that it would be interesting to look separately at configurations with fixed jmax instead of counting all 
eigenvalues belonging to an interval coming from all configurations with ja < jmax- 

Additionally it would be good to know the effect that our cutoff jmax = 50 has on the calculated eigenvalue 
spectrum. This is done in what follows. 

4.2.4 Lower Bound for the Spectrum 

For this purpose we will split the matrices to be calculated 13 into the sets of matrices indexed by configurations 
with fixed j'4 = j = jmax- We then inspect the ordered spectra of every set and try to find some regularity in 
the spin-configurations ji, ji, jz, producing the eigenvalues. 

As we have calculated all configurations up to jmax = 50 we have 100 sets of matrices Sj where each set is 
labelled by j = ... , jmax - 5, jmax- 

In each set we consider the first 100 positive eigenvalues ordered by values. Additionally to every eigenvalue 
we denote the spin configuration of the matrix giving rise to it and the position k the eigenvalue takes in the 
ordered list of eigenvalues of this set Sj. The thus achieved datasets are written into a file. 

Hence we create a function: 

A: (j,k)~ ♦A fc (j) (96) 

12 This maximal value is limited by to the capacities of the mathematical software Maple 7 on the computer used (Intel XEON 
machine with two 1.7 GHz processors). Future calculations will go much further, since Maple 7 is only an interpreter programming 
language. 

13 due to our setup 
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where A fe (j) is the fc-th eigenvalue in Sj. 

It turns out that the map (j,k) — > Xk{j) indeed displays a regularity, that is, the eigenvalues seem to 
produce series being separated from each other. Every series can be associated with a certain position in the 
ordered spectrum of a matrix set Sj with given y\ = j < jmax ■ The positions are taken from a minimal j on, 
as series k will not have contributions from Sj with too low j. It turns out, that the lowest eigenvalues of each 
matrix set Sj are precisely the lowest eigenvalues of the low dimensional matrices with low spin configurations 
ji < 32 < 33 < j < {ji + 32 + 32). 

Remarkably each of these matrices has rank smaller than 9 (that is, the nontrivial part of the characteristic 
polynomial can be reduced to a polynomial of degree less than or equal to 4), hence we can find analytic 
expressions for these lowest eigenvalues. 

We will give here a table containing the first 12 series of eigenvalues. In the second column we write down 
the smallest j from which (by inspection of the data) the noted order k is reached. Additionally we note the 
spin-configuration. The eigenvalues given are always the smallest ones/ of the according matrix giving rise 
to \k(j) with the given spin configuration. 

Surprisingly the first smallest eigenvalues are not equally distributed between even and odd configurations 
(the latter possessing O-eigenvalues) , but mainly contributed by the even configurations (we give in the second 
table the first odd configurations). Note again that these eigenvalues of Q v are the square of the eigenvalues of 
V v and that each eigenvalue has a multiplicity accordingly to the multiplicity table given above. 
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Even Configurations 



k 


valid j > 


3i 


J2 


33 


J4 






1 


i 

2 


i 

2 


1 
2 


3 


J 


Vi(i + 1) 


2 


2 


2 


1 

2 


1 


3-\ 




V0' + i)(2i-i) 


2^2 


3 


3 


3 
2 


3 
2 


3 




2\/ 17j 2 + 17? - 21 - \/208j 4 + 416j 3 - 344j 2 - 552j + 441 


2\/l7- ^208 


4 


13 
2 


1 

2 


3 
2 


3-1 


i 


2V3V0' + 1)0'-1) 


2^ 


5 


4 


1 

2 


2 


3-1 




2V2\/(i + l)(2j-3) 


4 


6 


7 
2 


1 


1 


J-l 




Vj(j + 1) 


4 


7 


20 


1 
2 


3 
2 


3 




2V(2j+3)(2j-1) 


4 


8 


n 

2 


3 
2 


2 


3~\ 




^ 108j 2 + 54j - 216 - Qy/228j 4 + 228j 3 - 903j 2 - 480j + 1152 


^108 - 6^228 


9 


65 
2 


5 
2 


5 
2 


3 


3 


too long but analytical expression 




10 


10 


1 
2 


5 
2 


J -2 


3 


2W0' + l)(?'-2) 


2^ 


11 


11 


1 

2 


3 


i-§ 


3 


2V3\/0' + l)(2j-5) 


2^6 


12 


18 


1 


3 
2 


i-f 


3 


2V3V0' + l)(2j-3) 


2^6 



Odd Configurations 



k 


ji 


J2 


33 


34. 


A* 00 


24 


1 


3 
2 


i-k 


j 


2^14i 2 + 7j - 16 


61 


2 


2 


j 


3 


2 y / 52j 2 + 52j - 114 - 18^4 j 4 + 8j 3 - 16j 2 - 20j + 33 



All of these expressions have to leading order the form 



lim **21 = Ck (97) 

j—>oa J 

where Ck increases with k. Thus for sufficiently large j the series j — > Xh(j) are approximate lines of different 
inclination. 

As an illustration we will give a plot of the first 8 eigenvalue-series, that is we plot for each j the first 8 
eigenvalues of the associated matrix-series Sj. Then we connect the first, second, third eigenvalues with a 
line. Here it becomes obvious, that for small j not all series are present, that is Xk(j) is ill defined below the 
threshold j given in the table for each k. 

Thus we are given a certain numerical criterion to decide, which part of the spectrum of the volume operator 
for the 4- vertex is already entirely calculated for a given cutoff jmax'- Given an eigenvalue A 2 , draw a horizontal 
line in figure [I)] and find the intersection with the first eigenvalue series that is, k = 1: 3 j with Ai(j) 2 = A 2 . 
The value j(A) at which this happens gives the maximal value jmax (A) which we have to consider in order to 
find configurations giving rise to eigenvalues < A, because all eigenvalues produced by j > j(X) are larger than 
A because numerically Xk{j) > Xi(j) V k > 1. 

According to the table above Ai(j) 2 = 2^/j(j + 1) = A 2 and therefore j(A) = — i + y§ + 4f • Thus for 

jmax = 50 we can trust to have computed the complete spectrum only for A < X ma x(jmax) = V% \J 'jmax (jmax + 1), 
i.e. X max = 1.4- V50 « 10. 
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Figure 5: The evolution of the first positive 8 eigenvalues Xk(j) ■— V 2 , k = 
1, .., 8 of Q v in dependence of j'4 = j. Note that each line represents eigenvalues 
with multiplicity given by the table in section T4. 2. II . 



4.2.5 Upper Bound for the Spectrum 

By observation of the numerical matrices in Sj it turns out, that the maximal eigenvalues \ max u) = (Vmlx) 2 
of configurations with fixed j'4 = j are contributed by matrices of the monochromatic vertex, that is ji = j'2 = 
js — j'4 = j as we expected from our estimates (I78H . The matrix elements of this special case we already wrote 
down in (JHZJ (0 < J12 < 2j): 



ra k (k = j'12) :=< ji2\qi23\ji2 - 1 > 



(2j + 1) 2 - (j 12 y 



(98) 



Now the theorem l|4.1|l provides us with upper bounds for the moduli of eigenvalues of a matrix in terms of 
its row- or column-sums. It is natural to ask now how the biggest eigenvalue Xmlix and the maximal row- or 
column-sum of the monochromatic matrix A of type l|70|l fit together. It is clear from the structure of A that 
the inequality given in theorem 14. II for the biggest eigenvalue reads: 



[|o fc | 



-I- |ofc + i| =: L 



k = l. 



Therefore we look for the maximal matrix element of A defined by 198H by differentiating the given expression 
with respect to j'12 and find the value ji2 ax of j'12 maximizing the matrix element. There are several extrema, 
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but the desired one turns out to be: 

1 



fi2 ax = ^V 24 -? 2 + 24 -?' + 12 + 6V'l6j 4 + 32j3+8j 2 -8j-2 



where for large j 



■ max J^oo 
3 12 — > 



a/3 



(99) 
(100) 



Since ji2 can only take positive integer values we then choose the maximal row sum which is given by 

Lmax = \dround(jY^ a:c ) I I ^round(j^ a:c ) — 1 1 

Plotting the quotient Xmax/L 

max cLS 3- function of — jmax 

we find, that this ratio converges to 1 as j increases: 




Vmax/Lmax 



10 20 30 40 50 



Figure 6: The quotient Xm ax /L max as a 
function of j 



Therefore we have numerical evidence for the following large j behavior of the biggest eigenvalue in a matrix 



set Si 



AmL J ^ L max {j) w 2\a } ™* (j5 aa! = j)! w \/3 y j 3 



(101) 



Here we have inserted jrf x in equation l|98^ for the matrix elements to obtain flj~' and approximated the 
maximal row sum L max by 2 ■ \ajrnax \. Finally we keep only the leading order of j in the expression and arrive 
at the result l|101[l . This coincides with the result obtained in |24j . 
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4.2.6 Spectral Density 



We now turn to a first investigation of how to get a reliable numerical estimate for the spectral density. 

Let us first briefly discuss the behaviour of the O-eigenvalues. We have proven that they only occur as a single 
member of the spectra of matrices with odd dimension. Therefore counting the number of odd configurations 
is equal to counting the number of O-eigenvalues. Since the total number of configurations grows with j^ al we 
also fit the total number of O-eigenvalues with a fourth order polynomial, whose coefficient in the leading order 
should approximately be half the coefficient we found when we fitted the total number of configurations, since 
we expect odd and even configurations to be nearly equally distributed (under restriction of 1881 . ). 



2e+07 - 



1.5e+07 - 



#0's 



1e+07 - 



5e+06 - 




Figure 7: The total number of O-eigenvalues con- 
tained in all configurations allowed by (|88J) de- 
pending on j max 



The fitted polynomial is obtained to be: 

#0 - eigenvalues(j ma x) = 3.33 • j ma x + 10.00 • j max + 10.' 



■ 2 
Jmax 



7.10 • jmax 



.85 



Indeed the coefficient 3.33 of j max is half the coefficient obtained for the total number of configurations before. 
The difference of the other coefficients seems to be caused by the restrictions given in 18811 . In what follows we 
will omit the number of 0-eigenvalues, since their behaviour does not contribute to the spectrum of non-zero- 
eigenvalues. Notice that their relative number as compared to the total number of all all eigenvalues is of the 
order of 



-i 

Jmax • 



Let us first recall how we define an eigenvalue density. We will take the square roots of the eigenvalues > of 
Q v obtained in the numerical computation and split the real axis labelling the eigenvalues V of V into identical 
intervals of the length AV. Then each eigenvalue V = n - AV (n = 1, 2, . . .) is assigned to an interval-number /„ 
defined by I n — round{^y + 0.5^ . In the third step we add up all eigenvalues belonging to the same interval 
/„ to get the number of eigenvalues in the interval [V — AV, V] . 



Now we define the Interval density Ni of eigenvalues in the interval / for fixed j\ = j by: 

#eigenvalues(I) 



NiU) ■= 



AV 



(102) 
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Since we want to have a normalized density p: 



pi dV 



we divide the interval densities Ni by the number of eigenvalues different from to get the final definition of 
the eigenvalue density: 

/•s ._ #eigenvalues(I) 

AV ■ (fetotal eigenvalues — #0 eigenvalues) 
These densities are then represented by a point at V = n ■ AV for each interval I n . These points are joined 
then (they can be fitted by polynomials for instance) to display the desired normalized eigenvalue density p. 
This gives us for instance for jn = j = 50 the following plot for a fit with a fourth order polynomial in the 
eigenvalues V (AV = 0.5): 



0.0035 



0.003 



0.0025 - 



0.002 - 



0.0015 - 



0.001 



0.0005 




Figure 8: The eigenvalue density for j'4 = 50 
(points) fitted by a 4th order polynomial (solid 
line). 



This can be done for every matrix set Sj with fixed j — j'4. Remarkably it seems to be true, that the 
eigenvalue densities in figEJare fitted quite well by 4th-order polynomials. But it is even more surprising that 
if we rescale the obtained densities, by putting their width W := Vmax(j) — Vmin{j) — > 1 and their height 
H(j) := max(fii(J)) — > 1 (where max(p/(j)) is taken from the fit curves) and plotting the resulting rescaled 
distributions for different values of 34 = j into the interval [0, 1] as given in figure lTTl and figure lT^l the distribution 
seem to possess a similar shape. 

Hence, the normalized distrbutions seem to be independent of j, that is, universal. This discussion suggests 
to try to define a limit distribution. By taking into account the behaviour of the ratio of the distance between 
the begin of the distributions V m in(j) to the value V(H(j)) at which the maximum H(j) of the distribution is 
situated and the total length of the distribution V m ax{j) — V m in{j). 

A(j) := YVm) - VminU) 

The ratio A(j) should tend to a constant value for j — j'4 — > 00 in the presence of a limit common shape of 
all distributions. Moreover we want to find out the quality of the fits taken by calculating the average squared 
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Figure 9: The eigenvalue densities for j'4 = 
10, 15, . . . , 50 (just look at the biggest eigenval- 
ues in order to identify, which curve belongs to 
which j). The points representing the eigenvalue 
density, are joined by lines. 



Figure 10: A fit of the spectra described in 
figure [5] by polynomials of 4th order. 
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difference between the fitted curves and the real spectra, given by (all quantities at given j): 



X 



1 



HVmam) 



m.SX.{pi)(y max — Vmin) 



I (fitted) 

2^ \PI- Pi 



i(v min ) 

These quantities seem to behave in a way which is convenient for us (see figures [T3l If 4H : 



0.22^ 

10 20 



]4=j 



Figure 13: The ratio A(j) dependent on ji = j 



Chi»2 



30 
j4=j 



Figure 14: The sum of the squared distances be- 
tween points of the eigenvalue densities and the 
fitted curves, % 2 , (defined as above) in depen- 
dence of J4 = j 



That is, the fit quality improves as j is growing and it seems possible, that the ratio A(j) indeed has a 
certain limit of ~ 0.25 . 



4.2.7 Density of the Eigenvalues of the Volume Operator - Discussion 

Despite the interesting properties and the occurrence of some systematics in the spectral density in the matrix 
sets, we were unable up to this point to give an estimate of the density of eigenvalues of the whole volume 
operator so far. 

The problem is first, that the 4-th order polynomials were chosen only for the reason, that they contain 
the lowest number of parameters (5), the achieved spectra can be satisfactorily fitted by. There has more 
computational work to be done to ensure the presence of la limiting eigenvalue distribution. If this turned out 
to be the case, then we have to look for certain points, the parameters of the fitted curves of the eigenvalue 
distribution are fixed by. 

What we know (at least our computations up to now encourage us to think that we know) are 4 parameters: 

3 i 

The maximum volume Vmax ~ jmax, the minimal volume V m i„ ~ jmax, the maximum of the distribution, 
situated (see fig. II 31 at ~ 0.25 • Vmax, and we know the total number E^ i= ^ ma, °' of eigenvalues , which can be 
obtained from the total number of eigenvalues E^ tot ^ for each j m ax by considering E^ i=:imax ^ = E' tot > (j max ) — 
E^ tot ' {jmax — |). Then we could compute the eigenvalue density in an interval [Vi, V2], by a superposition of 
all configurations with Vmin < V2 and Vmax > Vl. 

The second problem is then to give a fit-function for the total density of eigenvalues for all configurations. 
There is one conjecture, based on a similar behaviour for the spectrum of the area operator (which is easier to 
handle), that this eigenvalue density should behave as 

p{V) = 

We have tried to fit the part of the spectrum we have fully calculated (V < V m ax{jmax = 50)) by this formula. 
But it seemed to be impossible to give certain values to the three parameters, especially to 7. Maybe the 
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calculated part of the spectrum is still to small, i.e. not sufficient for statistics or one must consider also higher 
valent vertices. Of course, the conjecture could also be wrong. So it is left as an open (but inspired by our done 
calculations not hopeless) task, to fix the density of eigenvalues. 

Let us conclude by displaying here the complete calculated part of the eigenvalue density of the volume oper- 
ator for jmax = 50, according to our numerical criterion, that we can rely on for given j ma x on the part of the 

spectrum with V < J 2y/ j m ax(jmax + 1). That is for jmax = 50 (as we calculated) we have obtained the full 
spectrum up to V ~ 10. 

By inspection of figure we can extend this part up to V ~ \/200 ~ 14 if we draw a horizontal line into 
figure at a given 10 < V < 14 and count the eigenvalue series situated below that line. If we assume, that 
these series grow linearly with growing j m ax, as we expect, and there do not occur additional eigenvalue series at 
higher j m ax then we can simply extend the curves j — > Afc(j) linearly for k = 1 . . . 8 and can thereby estimate 
(approximately) the additional contributions not yet calculated explicitly for V 2 < 200 that is for the k = 7 
eigenvalue series. 

Therefore we display the original spectrum (circles) and the extended spectrum (boxes) normalized with 
respect to the total number of non-zero eigenvalues of the original spectrum (we have chosen again an interval 
width AV = 0.5, the eigenvalue density is defined as in 1)103^ ) in figure IT51 
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Figure 15: The eigenvalue density of the total volume operator (circles) and the 
extended eigenvalue density (boxes) for j max = 50 on a 4-valent vertex 
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5 Summary &; Outlook 



In this article we have analyzed the spectral properties of the Volume Operator defined in Loop Quantum 
Gravity. 

We discussed the matrix representation of the volume operator with respect to gauge invariant spin network 
functions and were able to derive a drastically simplified formula for the matrix elements of the volume operator 
with respect to the latter. 

It turned out that there exist certain selection rules for the matrix elements and all the matrices are i-times 
an antisymmetric matrix with the structure of a Jacobi matrix, meaning that non-vanishing matrix elements 
are only situated on certain off- diagonals. 

We were able to determine the kernel, that is, the eigenstates for the eigenvalue 0, of the volume operator 
with respect to the gauge invariant 4- vertex analytically as given in l|85|l . I|86|l . We have done numerical inves- 
tigations for the gauge invariant 4-vertex. Our numerical investigations support the analytical estimate, that 



there exists a smallest eigenvalue V m in dependent of the the maximal spin j max via Vmin ~ ( jmax J and a 



maximal eigenvalue V m ax ~ [Jmax J ■ Therefore we were able to find certain numerical indicators for the com- 
pleteness of a numerically computed part of the spectrum. Moreover, we found that the geometrical intuition 
is reflected in the spectrum: At given jmax the lowest non - zero eigenvalues come from rather "distorted", 
almost flat tetrahedra with large spin on some edges and low spin on the others. On the other hand, the largest 
eigenvalues come from regular tetrahedra with large spin on all edges. 

For future analysis one should extend the numerical calculations for the gauge invariant 4-vertex and higher n- 
valent vertices to verify and possibliy extend (for higher valence vertices) the regularities of the spectra obtained 
(for the gauge invariant 4-vertex) for single matrix sets Sj with j'4 = j at higher spins. In particular it would 
be interesing to see whether there exists a volume gap as n — > 00. This, however, requires more computing 
power and better programming than we have used in this paper. 

The formula derived for the matrix elements with respect to a gauge invariant n-valent vertex can be used 
to analyze the whole spectrum of the volume operator numerically and analytically. Further simplifications are 
conceivable. 

Notice again that the restriction to 4-valued vertices and j ma x < 50 was only due to the computational ca- 
pacity of the used mathematical software MAPLE 7 and computer. By using a compiler-based programming 
language and optimized numerical matrix-diagonalizing routines we expect to be able to go much beyond the 
computational limits above. 

Due to the results presented in this paper it seems to the authors, that there are good chances for getting 
sufficient control about the spectral behaviour of the volume operator in the future, especially when it comes 
to dynamical questions in LQG. 

As a first qualitative application in that respect, notice the following: 

We have shown analytically that the volume operator of full LQG has zero eigenvalues at arbitrarily large jmax 
and that their number grows as j max as compared to the total number of eigenvalues which grows as j max , at 
least for the gauge invariant four vertex which should be the most interesting case from a triangulation point 

1/2 n 

of view. Moreover, the volume gap increases as jmax- It follows that the full spectrum contains many "flat 
directions" or "valleys" of zero volume and the walls of the valleys presumably get steeper as we increase jmax ■ 
Therefore we might find arbitrarily large eigenvalues as close as we want to zero eigenvalues and hence the 
"derivative" (rather: difference) of the spectrum around zero volume which enters the Hamiltonian constraint 
through the curvature operator, while well - defined as shown in 4 5 , could be unbounded from above. There- 
fore, the full spectrum could not share an important property of the spectrum in the cosmological truncation 
of LQG 1 1 12 1 where the derivative of the spectrum at zero volume is bounded from above. As this property has 
been somewhat important in |12|. some of the results of |12| might have to be revisited in the full theory. The 
challenge would be to show that the curvature expectation value remains bounded when the system is prepared 
in a semiclassical state for LQG, see e.g. |M5| and references therein. First evidence for this and further analysis 
will be presented soon in I36| thus reaffirming the spectacular results of |12| in the full theory. 
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A Basics of Recoupling Theory 

A.l Angular Momentum in Quantum Mechanics 

For the angular momentum operator J = (Ji, J2, J3) (where each component has to be seen as an operator) we 
have the following commutation relations: 

Ui , Jj] = i ■ djk Jk [j 2 , Jj] = (104) 

with j 2 = j 2 = Jl + jf + Ji. 

Additionally we can define 

J+ = Ji + iJ 2 J- = Ji - iJ 2 (105) 

with (using l|l()4fl ): 

[j 2 , J±] = [j 3 ,J+] = J+ [i 2 , J-] = ~J- [j+,J-] = 2J 3 (106) 

Since every angular momentum state is completely determined 14 by its total angular momentum quantum 
number j (where J 2 — j(j + 1)) and one component say J3 (where J3 = — j, —j + 1, ... ,j — l,j ) we then 
associate for certain j a 2j + 1 dimensional 15 Hilbert space TL equipped with an orthonormal basis \j m >, 
m = Ji where: 

< j m\j m >= 5 mm > (107) 

The I j m > simultaneously diagonalize the 2 operators of the squared angular momentum J 2 and the 
magnetic quantum number J3 |14|: 

J 2 \j m >= j(j + m > Ja\j m >= m\j m > (108) 

That is, \ j m > is a maximal set of simultanous eigenvectors of J 2 and J3. 
On these eigenvectors the other operators act as 

J+\j m >= V J(j + 1) - m(m + 1) |j m + 1 > J-\j m>= + 1) - m(m - 1) | j m - 1 > (109) 



A. 2 Fundamental Recoupling 

Equipped with a small part of representation theory we can easily understand what happens if we couple several 
angular momenta. For that we first repeat the well known theorem of Clebsh & Gordan on tensor products of 
representations of SU(2): 

Theorem A.l Clebsh & Gordan 

Having two irreducible representations Hj 1 , 7r J2 of SU (2) with weights j\ , j 2 their tensor product space splits into 
a direct sum of irreducible representations Hj 12 with \ j\ — ja| < J12 < ji + ji such that 



Equivalently we can write for the resulting representation space H}- D ^ = li}- Dl ' ® li} D2 ' (where D\ = 
2ji + 1, D2 = 2j2 + 1, D — Di ■ D2 denote the dimensions of the Hilbert spaces): 

h {d) =h (Di) ®h (D2) = 2 H (2jl2+1) (110) 

J12 = I J l —31 

Or in other words: If we couple two angular momenta we can get resulting angular momenta j'12 

varying in the range |ji - j 2 \ < J12 < ji + h- 

The tensor product space of two representations of SU(2) decomposes into a direct sum of representation 
spaces, that is one space for every possible value of recoupling j'12 with the according dimension 2ji2 + 1. 

In the sense that we have a maximal set of simulaneously measureable observablcs. 
for fixed j there are 2j + 1 values which J3 can take. 
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A. 3 Recoupling of Two Angular Momenta 

According to lllUl we can expand for each value of 312 the elements \j\ {31,32) M > G H (2jl2+1) called 

'coupled states' into the tensor basis \ji mi > ® |j 2 m 2 > of T6 D ^: 

\h 32\ 3Mji, 32), M > = < ji mi;j 2 m 2 \ji 32; 312(31, 32), M > \j\ mi > ®\j 2 m 2 > (111) 

mi +m 2 — M v 

°™»l m 2 

Here C mim2 G 1R denotes the expansion coefficients, the so called Clebsh-Gordan coefficients. On the right hand 
side \ji mi; 32 m 2 >= \ji mi > <g> jj 2 vn,2 > 

If we change the order of coupling then C mim2 changes its sign: 

< ji mv,h "i2|ji2(ji, 32) M = m 1 +m 2 > = (-i)^2-Ji-J2 < j 2 m 2 ;ji m 1 \j 12 (j2, ji) M = m 1 + m 2 > 

= (-1)-Ji2+Ji+J2 < J2 m2 .j 1 mi|ji 2 (j 2 , Ji) M = mi + m 2 > 

As =Fii2 ± (ii + Ja) is an integer number we are allowed to switch the signs in the exponent of the factor (— 1). 
The coupled states again form an orthonormal basis: 

< 3i 32; 312(31, 32), M\n j 2 ;ji2{ji,j2), M > = 1 (112) 

< ji 32 ; 312 (ji ,32), M\ji 32 ; 312 (ji ,32), M > = S~ ijia 8% M (113) 



In piHI 5~ 2 j comes from the orthogonality of the Hy 2 i'i' 2+1 ) i n l|f f 5jj m is caused by the othogonality 
of the single \ji mi >, \ j% m 2 > l|107|l - since always M = mi + m 2 . 

Normalization of the recoupled states 11121 implies according to 110711 : 

I < k h; 312(31, 32), M\j! 32; 312(31, h), m > I 2 = l^i^l 2 = 1 (114) 

m l+ m 2~ M m l+ m 2~ M 

Furthermore the Clebsh — Gordan — coef ficients are all real, which is not obvious, but a result of two 
conventions one usually requires |14| : 

1) \ji 32; 312(31, h) = ji +32 M = ji + 32 >= \ji mi = ji > ® |j 2 m 2 = j 2 > 

2) All matrix elements of J^ Dl \ which are nondiagonal in \ji 32', 312(31, 32) M > are real and nonnegative. 

The maximal set of simultanuosly diagonalizeable (that is commuting) 2 ■ 2 operators l|108> of the single 
Hilbertspaces Tl^ Dl \7i^ D2 ^ is then in JiP 3 ^ replaced 16 by 4 operators: total angular momentum (J^ D ^) 2 , total 
projection quantum number J^ D \ single total angular momenta (J^ Dl ^ , (J^ D2 ^) 2 : 

{J (D }) 2 = {j(Dl) + J( D 2)) 2 j(D) = j( Dl) + j( D 2) {j(Dl)) 2 {j( D 2)) 2 (ng) 

which are simultanously digonal in the new basis manifested through: 

(J {D) f \ji 32; 312(31, 32), M > = ji 2 0'i2 + l) \ji 32; 312(31, 32), M > 

4 D) \ji hi 312(31, 32), M > = M \ji 32; 312(31, 32), M > 

(J (Dl) ) 2 \ji 32; 312(31, 32), M > = ji(ji + 1) \ji 32; 312(31, 32), M > 

(J [D2) ) 2 \ji 32; 312(31, 32), M > = 32(32 + 1) \ji 32; 312(31, 32), M > (116) 

A. 4 Recoupling of Three Angular Momenta - 6j-Symbols 

In this way we can expand the recoupling of three angular momenta in terms of CGC. 

\3i2(ji, h), j(ji2, h) > = \ji 2 (ji,j2) js; 3(312, js.) M = mi +m 2 + m 3 > 

= Y, < • ?12 m i2',h m 3 \ji2 33; 3 mi2 + m 3 > • 

m 12 m 3 

\ji2 mi2 > |ja rns > 



16 We are a bit sloppy in using this notation, correctly we would have to write: 
(J(0))2 = (jCDi) g, lw(D2) + lw(Di) J(D 2)) 2 j(D) = j(D t ) ^ 1 ^ (D2) + 1 ^ (Di> ^ j(D 2 ) ^ 

{ j( Dl)) 2 = (J (D l))2 ^ lw(jJa) (j(C2 ) ) 2 = 1h(Di)0(j (D 2)) 2 
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= ^2 < 312 m 12 ;j3 m 3 |ji2 jy,j m 12 + m 3 > • 

"112 "13 

22 < h m l>h m 2\ji h;ji2 TO12 = mi + m 2 > ■ 

mi 771-2 

|ji mi > \j 2 m 2 > |j3 m 3 > 
= 2~2 < to mi + rn.2; js M — mi — rna\ji2 ja; j M > ■ 

rn\ 7Tt2 

< ji mi; 32 m 2 \ji 32; 312 mi + m 2 > ■ 

\ji mi > \j 2 m 2 > |j3 M - mi - tri2 > (117) 

As we can see, as we couple angular momenta successively, the order of coupling plays an important role. 
Different orders of coupling will lead to different phases of the wavefunctions (see 111211 ). Concerning this it 
would be nice to have a transformation connecting different ways of recoupling. This tranformation between 
two different ways of coupling 3 angular momenta 31,32,33 to a resulting j defines the 6j-symbols, see appendix 

m 

A. 5 Recoupling of n Angular Momenta - 3ny-Symbols 

As mentioned before the case of successive coupling of 3 angular momenta to a resulting j can be generalized. 
For this purpose let us first comment on the generalization principle before we go into detailed definitions. 

Theorem rO can be applied to a tensor product of n representations ■Kj 1 ® irj 2 ® . . . ® 7Tj„ by reducing out 
step by step every pair of representations. This procedure has to be carried out until all tensor products are 
reduced out. One then ends up with a direct sum of representations each of them having a weight corresponding 
to an allowed value of the total angular momentum the n single angular momenta 31,32, ■ ■ ■ ,jn can couple to. 

But there is an arbitrariness in how one couples the n angular momenta together, that is, the order by which 
7Tjj ® 7Tj 2 ® ■ ■ ■ ® TTj n is reduced out (by applving lA.lfl matters. 

Let us now have a system of n angular momenta. First we fix a labelling of these momenta, such that we 
have 31,32, ■ ■ ■ ,jn- Again the first choice would be a tensor basis \j rh > of all single angular momentum states 
\jk rrik >, k = 1 . . . n defined by: 

n 

\j rh >= \(ji,j 2 , ■ ■ ■ ,jn) (mi,m 3| .. . ,m„) >:= \j k m k > (118) 

k-l 

with the maximal set of 2n commuting operators (Ji) 2 , Jf, (I = 1, . . . , n). 

Now we proceed as in section TA.3I finding the commuting operators according to 1116H . that is a basis in 
which the total angular momentum (Jtot) 2 = (J) 2 = {Ji + J 2 + . . . + J n ) 2 is diagonal (quantum number j) 
together witch the total magnetic quantum number jf ot = J 3 = jf + J 2 + . . . + J„ (quantum number M). 

As (J) 2 and J 3 are 2 operators, we need 2(n — 1) more quantum numbers of operators commuting with 
each other and with (J) 2 and J 3 to have again a maximal set. We choose therefore the n operators (Ji) 2 , 
I — 1, ... ,n of total single angular momentum (quantum numbers (ji, . . . ,j n ) := j). So we are left with the 
task of finding addtional n — 2 operators commuting with the remaining ones. For this pupose we define: 

Definition A.l Recoupling Scheme 

A recoupling scheme \g(IJ) j j m > is an orthonormal basis, diagonalizing besides (J) 2 , J 3 , (Ji) 2 
(I = 1, . . . , n) the squares of the additional n — 2 operators G 2 , G3, . . . , G n -i defined as 17 : 

Gi := Ji, G 2 := Gi + Jj, G3 := G 2 + Ji, Ga '■= G3 + J 2 , Gi := G1-1 + Ji- 2 , 

G1+1 := Gi + J1-1, G1+2 := Gj+i + J1+1, G1+3 := Gi+ 2 + J1+2, Gj := Gj-i + Jj-i, 

Gj+l := Gj + Jj+l, Gj+2 '■= Gj+l + Jj + 2, • • • , G n -1 := G n - 2 + Jn-l 

The vector g(IJ) := [g 2 {ji , jj), 53(32, ji), ■ ■ ■ , 31+1(9/ , ji-i), gi+2{gi+i, ji+i), ■ ■ ■ , g.j(gj-i,jj-i), 
i(gn-2, jn-ifj carries as quantum numbers the n — 2 eigenvalues of the operators 

(G 2 ) 2 , (G„_i) 2 

So we recouple first the angular momenta labelled by 7, J where I < J and secondly all the other angular 
momenta successively (all labels are with respect to the a fixed label set), by taking into account the allowed 

Note, that formally G n ■= G n -± + J n = Jtotal 
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values for each recoupling according to theorem I A. II 



Let us define furthermore the so called standard recoupling scheme or standard basis: 
Definition A. 2 Standard Basis 

A recoupling scheme based on the pair (I, J) — (1, 2) with 

K 

Gk = Jl 

L=l 

is called standard basis. 

Using definition IA.1I with the commutation relations 11U4H and the fact, that single angular momentum 
operators acting on different single angular momentum Hilbert spaces commute 18 , one can easily check that for 
every recoupling scheme 

(i) the G/'s fulfill the angular momentum algebra l|104[l . 

(ii) (J) 2 , (J/) 2 , (G K ) 2 , J 3 commute with each other V /, K = 1 ... n 

Note, that it is sufficient to prove these two points in the Standard basis £(12), because every other basis g(U) 
is related to it by simply relabelling the n angular momenta. 

We have thus succeeded in giving an alternative description of an n angular momenta system by all possible 
occurring intermediate recoupling stages Gi instead of using the individual magnetic quantum numbers. 

Obviously every orthonormal basis spanned by a recoupling scheme \g(IJ) j j m > is singled out by the 
labelling, namely the index pair (7J) and therefore not identical, as we have already seen in the case of the 
two angular momentum problem. So we are in need of a transformation connecting the different bases that is 
expressing one basis, e.g. belonging to the pair (IJ), in terms of another basis, e.g. belonging to the pair (KL), 
respectively. This leads to the following. 

Definition A. 3 3nj -Symbol 

The generalized expansion coefficients of a recoupling scheme in terms of the standard recoupling scheme are 
called 3nj -symbols: 

| g(IJ) j jm>= <g (12) Jjm \ g(IJ) jjm> | £'(12) j j m > 

The summation has to be extended over all possible values of the intermediate recouplings 
g'(12) — (fl^C/i, J2)) 33(52; J3)) •••) 9n-ii,9n-2i jn—l) t that is all values of each component g' k allowed by theo- 
rem TA.l\ 

In calculations we will supress the quantum numbers j,j,m, since they are identical all the time, and write 
abbreviating < g(U) \ £(12) >. Note additionally, the followoing properties of the 3nj'-symbols: 

(i) They are real, due to the possibility to express them as Clebsh-Gordan-coefficients: 

< g(IJ) I £(12) > = < £(12) I £(/J) > 

(ii) They are rotationally invariant, i.e. independent of the magnetic quantum numbers mu occurring in 11181 . 

That is J J j\ = whenever I ji J 
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B Properties of the 6j-Symbols 

In this section we will give an overview on the 6j-symbols because they are the basic structure we will use in 
our recoupling calculations, every coupling of n angular momenta can be decomposed into them. For further 
details we refer to |14| . |15) . 

B.l Definition 

The 6j-symbol is defined as |14|.d 92: 



h i £ } ; [{2jl2 + 1)(2j23 + i)]^(-i) J1+32+J3+J 

X < J12 (jl , 32) , j (jl2 , js ) b'23 O'a ,33), j (jl , J23 ) > 

= [(2j 12 + l)(2j23 + l)]-4(-l)«+» + »+' 

x ^2 < ji m 1 ;j 2 m 2 \ji 32 ji2 mi+ma> 

mi 7712 

x < ji2 ra x + m 2 ;jz m — m\ — m 2 |ji 2 jz j m> 
x < j 2 m 2 ; jz m — mi — rn 2 \j 2 jz 323 m — mi > 

x < ji m 1 ;j 2 z m - mi\ji j 23 j m> (119) 
The terms under the summation are called Clebsh-Gordon-CoefEcients. 

B.2 Explicit Evaluation of the 6j-Symbols 

A general formula for the numerical value of th3 6j-symbols has been derived by Racah [El, |H|, p. 99: 

31 32 312 1 = A(j 1 ,j 2 , Jl2 )A(il,j,j 2 3)A(i3,j2,j23)A0 3 ,J,,jl2v{f J ? f» 
J3 3 32Z I I J3 J J23 



A(o, 6,c) = 



(0 + 6 - c)!(a - 6 + c)!(-a + 6 + c)! 



(a + 6 + c+l)! 



x [(n - ji ~ ia - ji2)!(n - ji - j - j 2 3)!(n - jz - j 2 - j 2 z)K n ~ 33 - j - jw)!]" 1 x 

x [(ji + j 2 + jz + j - n)!(j 2 + j i2 + j + j 23 - n)!(ji2 + ji + j 23 + js - (120) 

The sum has to be extended over all positive integer values of n such that no factorial in the denominator has 
a negative argument. That is: 

max[ji + j 2 + ji2 , ji + j + j23 , jz + 32 +j2z , jz +j +312] < n < min[ji + j 2 + j 3 +j, j 2 + ji 2 + j + j 2 3 , ji2 + ji + j23 +jz] 

Remark From fT^Hi we are provided with some additional requirements, the arguments of the 6j-symbols have 
to fulfill: Certain sums or differences of them have to be integer for beeing proper(=integer) arguments 
for the factorials: 

from A(a,b, c) one gets: 

• a,b,c have to fulfill the triangle inequalities: (a + b — c) > 0, (a — b + c) > 0, (—a + b + c) > 0, 

• (±a ± b ± c) has to be an integer number 
from the w-coefncient one gets: 

• ji + 32 + jz + j, 32 + ji2 + j + j2z , ji2 + ji + j 2 3 + jz are integer numbers. 

The following (trivial but important) relations are frequently used in calculations involving 6j-symbols: 

{-if = (-iy z Vz e z 

(~lf z = 1 Vz 6 Z 

(_1)3* = Vfc=~ withzeZ (121) 
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B.3 Symmetry Properties 

The 6J-symbols are invariant 

• under any permutation of the columns: 
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• under interchange of the upper and lower arguments of two columns at the same time. E.. 
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B.4 Orthogonality and Sum Rules 

Orthogonality Relations 

y (2j la + l)(2j 12 + 1) ( jl il2 } ( '} •': 3 ' 2 \ =6. , 

^ J A ^ \ J3 J J23 )\j 3 3 323 J J12J 12 

Composition Relation 

y^C_ 1 y23+J31+J12 ( - 272 + 1) J Jl 32 312 \ { 32 33 323 \ = { 33 31 331 \ 

^ \ 33 3 323 J \ ji 3 33i J \ 32 3 ji2 J 



Sum Rule of Elliot and Biedenharn 
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xV(-ir (2ji24 + 1) I i3 •': j2A ] { ■'• jl j . 12 } i jl2 \ 
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C Comment on the Smallest Non- Vanishing Eigenvalue 



In this section we will briefly summarize, what can be done to obtain a lower bound of the spectrum of the 
matrices occurring when expressing the Volume Operator on a recoupling scheme basis at the gauge invariant 
4- vertex. This is mainly done to illustrate the remarkable symmetries in that case. The idea is to obtain a lower 
bound of the non-zero-eigenvalues by applying theorem 14.11 . on the inverse matrix, giving an upper bound for 
its eigenvalues and therefore a lower bound for the non-zero volume spectrum. 

The general form for the gauge-invariant 4-vertex was obtained in 1701 

/ -ai \ 

<2l U — 0,2 '■ 



A = 





a 2 







V o 



d-n-l 







(127) 



We have explicitly discussed the 0-eigenvalues contained in the spectrum of A in section 14.1.21 We know, 
that only in the odd-dimensional case the matrix A is singular containing one 0-eigenvalue with the according 
eigenvector 



(n). 
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(128) 

0204 ■ . . . • a„_i 

where n := dim A and x an arbitrary scaling factor. We will denote the k th element of ' n *S! by ' n '9k- F° r 
technical reasons we will set x = r „ in the following to obtain: 



S :=y 
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(129) 



C.l Even Dimension n of A 

Since A is regular in that case we can invert it to find: 
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where we used the components 
with scaling factor y = 1. 
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(130) 



Hfc of the O-eigenvector X S of the (n — l)-odd dimensional case 
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C.2 Odd Dimension of A 

Since A is not regular we have to project out its nullspace ' n fe first (with arbitrary prefactor y). That can be 
done by applying a similarity transformation W on A to obtain R := W~ 1 AW: 
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■ 1 d n-1 
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Now one can check, that (all even components of vanish, again y = 1): 



(131) 
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(132) 
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To discuss the non-zero-spectrum of A it is sufficient to discuss the now regular submatrix R being the submatrix 
of R with n th row and column deleted. One obatins: 
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(133) 
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Here Mij = 



and R(ij) is a shortcut for the submatrix of R one obtains by deleting row i and column 
j. Basically this is the definition for the matrix element of the inverse matrix. Since unfortunately the Mij are 
hard to control (but of order -) we are unable to give an explicit upper bound for the spectrum of RT 1 and 
therefore a lower bound on the spectrum of A according to theorem l|4.1|l . It is remarkable, however, that half 
of the sructure of i jlHOI is being reproduced by the odd dimensional case. 



54 



References 

[1] C. Rovelli, L. Smolin, "Discreteness of volume and area in quantum gravity" Nucl. Phys. B 442 (1995) 
593, Erratum : Nucl. Phys. B 456 (1995) 734 

[2] A. Ashtekar, J. Lewandowski 

"Quantum Theory of Gravity II: Volume Operators", 

[arXiv: |gr-qc/971103T] 
[3] C. Rovelli, 

"Quantum Gravity" , Cambridge University Press 2004, at press. 
"Loop Quantum Gravity", Living Rev. Rel. 1 (1998) 1, 
[gr-qc/9710008| 
T. Thiemann, 

"Modern Canonical Quantum General Relativiy", 
to appear in Cambridge University Press, 
[arXiv: |gr-qc/01 10034] . 

"Lectures on Loop Quantum Gravity", Lecture Notes in Physics, 631 (2003) 41 - 135, 
gr-qc/0210094 

A. Ashtekar, J. Lewandowski, 

"Background Independent Quantum Gravity: A Status Report" , 
gr-qc/0404018 
[4] T. Thiemann, 

"Anomaly-free Formulation of non-perturbative, four-dimensional Lorentzian Quantum Gravity" , Physics 
Letters B380 (1996) 257-264. 
gr-qc/9606088 

[5] T. Thiemann, 

"Quantum Spin Dynamics (QSD)", Class. Quantum Grav. 15 (1998) 839-73, gr-qc/9606089 ; "II. 
The Kernel of the Whcclcr-DcWitt Constraint Operator", Class. Quantum Grav. 15 (1998) 875-905, 
gr-qc/9606090 ; "III. Quantum Constraint Algebra and Physical Scalar Product in Quantum General 
Relativity", Class. Quantum Grav. 15 (1998) 1207-1247, |gr-qc/9705017] ; "IV. 2+1 Euclidean Quantum 
Gravity as a model to test 3+1 Lorentzian Quantum Gravity", Class. Quantum Grav. 15 (1998) 1249-1280, 
gr-qc/9705018 ; "V. Quantum Gravity as the Natural Regulator of the Hamiltonian Constraint of Matter 
Quantum Field Theories", Class. Quantum Grav. 15 (1998) 1281-1314, [gr-qc/9705019 ; "VI. Quantum 
Poincare Algebra and a Quantum Positivity of Energy Theorem for Canonical Quantum Gravity", Class. 
Quantum Grav. 15 (1998) 1463-1485, gr-qc/9705020 

[6] T. Thiemann, 

"The Phoenix Project: Master Constraint Programme for Loop Quantum Gravity", 
gr-qc/0405080 

[7] R. Loll, Phys. Rev. Lett. 75 (1995) 3048 

[8] R. De Pietri, 

"Spin Networks and Recoupling in Loop Quantum Gravity", Nucl. Phys. Proc. Suppl. 57 (1997) 251 
[arXiv: |gr-qc/970104l1 

[9] R. De Pietri, C. Rovelli, 

"Geometry Eigenvalues and Scalar Product from Recoupling Theory in Loop Quantum Gravity" , Phys. Rev. 
D54 (1996) 2 664 
[arXiv: gr-qc/9602023 v2] 

[10] T. Thiemann, 

"Closed formula for the matrix elements of the volume operator in canonical quantum gravity", J. Math. 
Phys. 39 (1998) 3347-3 371 
[arXiv: gr-qc/9606091 . 
[11] J. Brunnemann, 

"Spectral Analysis of the Volume Operator in Canonical Quantum General Relativity", Diploma Thesis, 
Humboldt Universitat zu Berlin, December 2002 

[12] M. Bojowald, H. A. Morales - Tecotl, 

"Cosmological Applications of Loop Quantum Gravity" , 

gr-qc/0306008 
[13] A. Perez, 

"Spin Foam Models for Quantum Gravity", Class. Quant. Grav. 20 (2003) R43 
gr-qc/0301113 



55 



[14] A. R. Edmonds 

"Angular Momentum in Quantum Mechanics," Princeton University Press, Fourth printing 1996 

[15] D. A. Varshalovich, A. N. Moskalev, V. K. Khersonskii, 

"Quantum Theory of Angular Momentum" , World Scientific, 1988 

[16] R.U. Sexl H.K. Urbandtke, 

"Relativitat, Gruppen, Teilchen", Springer Verlag Wien/New York 1982, Zweite, erweiterte Auflage 

[17] S.M. Carroll, 

"Lecture Notes on General Relativity", 
[arXiv: lgr-qc/9712019] 

[18] Robert M. Wald, 

"General Relativity", The University of Chicago Press, 1984 

[19] G. Racah, "Theory of Complex Spectra. II", 

Phys. Rev.62 (1942), published in [D] 
[20] L.C. Biedenharn H. van Dam, 

"Quantum Theory of Angular Momentum" (a collection of reprints and original papers edited by the 

authors), Academic Press New York and London 1965 
[21] Alfred Horn 

"Eigenvalues of sums of Hermitian matrices", Pacific J. Math. 12 (1962), 225-241 

[22] William Fulton 

"Eigenvalues, invariant factors, highest weights and Schubert calculus", 
[arXiv: math.AG/9908012| 

[23] D.E. Holz H. Orland A. Zee, 

"On the Remarkable Spectrum of a Non-Hermitean Random Matrix Model" , 
arXiv:math-ph/0204015 

[24] M. Seifert, 

"Angle and Volume Studies in Quantized Space", 

[arXiv: gr-qc/0108047 
[25] F.J. Dyson, 

"The Dynamics of a Disordered Linear Chain", Phys. Rev. 92, 1331 (1953) 

[26] F. R. Gantmacher 

"Matrizentheorie" , VEB Deutscher Verlag der Wissenschaften, Berlin 1986 

[27] R. Zurmuhl, 

"Matrizen", Springer Verlag Berlin/Gottingen/Heidelberg 1964, 4. Auflage 
[28] M. Marcus H. Mine, 

"A Survey of Matrix Theory and Matrix Inequalities", Dover Publications, Inc. New York 1992 

[29] H. Grassmann 

"Algebra und Geometrie - Vorlesungsscript" 

http: / /www. irm-mathematik.hu-berlin.de/~hgrass 

[30] P.A.M. Dirac, "Lectures on Quantum Mechanics", 

Academic Press INC. New York/London, Second printing 1967 

[31] W. Nolting, 

"Grundkurs Theoretische Physik: 5 Quantenmechanik, Teil 2: Methoden und Anwendungen" , Friedrich 
Vieweg & Sohn Verlgsgesellschaft mbH, Braunschweig/Wiesbaden, 1997 

[32] N.N. Bogoljubov, D.V. Sirkov, 

"Quantenfelder" , Physikverlag Weinheim / VEB Deutscher Verlag der Wissenschaften, Berlin 1984 
[33] A. Wintner, 

"Spektraltheorie der unendlichen Matrizen", Verlag von S. Hirzel in Leipzig 1929 
[34] Teubner-Taschenbuch der Mathematik, Bd. 1+2, B.G. Teubner Stuttgart, Leipzig 1996 
[35] A. Ashtekar, J. Lewandowski 

"Relation between Polymer and Fock Excitations", Class. Quant. Grav. 18 (2001) L117-L128 

gr-qc/0107043 

T. Thiemann 

"Complexifier Coherent States for Quantum General Relativity", gr-qc/0206037 

[36] J. Brunnemann, T. Thiemann, 

"Towards the Cosmological Sector of Loop Quantum Gravity" , in preparation 



5G 



